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ABSTRACT 


It  is  difficult  to  propagate  a  diffraction-limited  laser 

beam  through  the  atmosphere,  since  the  atmosphere  contains 

random  index  of  refraction  fluctuations.  Two  parameters  that 

characterize  the  atmosphere  for  optical  propagation  are  the 

atmospheric  isoplanatic  angle,  0  ,  and  the  refractive 

2 

turbulence  structure  parameter,  Cn  .  This  dissertation  deals 

2 

with  improved  methods  for  measuring  9Q  and  Cn  profiles  using 
optical  techniques. 

By  apodizing  the  receiver  telescope  aperture,  one  can 

improve  the  weighting  function  for  isoplanatic  angle 

measurement  substantially  over  previous  systems.  We  find 

that  the  weighting  function  is  not  significantly  affected  by 

inner  scale  changes  with  altitude  and  that  the  error  in 

isoplanatic  angle  measurement  from  strong  low  altitude 

turbulence  (z  <  1  km)  with  this  weighting  function  is  small. 

Data  collected  with  the  improved  isoplanometer  shows  temporal 

trends  in  the  isoplanatic  angle  on  the  order  of  90  seconds 

that  have  not  been  observed  before. 

Direct  inversion  of  the  amplitude  covariance  function 

(including  aperture  averaging  effects)  to  yield  refractive 

turbulence  profiles  is  known  to  be  ill-posed.  I  suppress 

this  condition  using  Tikhonov  regularization  and  reproduce 

2 

refractive  turbulence  profiles  from  actual  C  data  with  some 

"  £3 


success . 
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I.  INTRODUCTION 


If  one  wishes  to  propagate  a  spatially  and  temporally 
coherent  electromagnetic  wave  through  the  atmosphere  or  any 
random  medium  without  distortion,  one  must  somehow 
compensate  for  the  random  phase  and  amplitude  perturbations 
induced  by  the  medium.  Adaptive  optics  is  the  field  of 
physics  that  attempts  to  compensate  for  propagation  through 
the  atmosphere,  principally  by  deforming  mirrors  and 
unconventionally  by  using  nonlinear  optical  materials.  In 
order  to  successfully  conjugate  the  effects  of  the 
atmosphere,  knowledge  of  certain  atmospheric  parameters  is 
crucial.  These  parameters  are  the  isoplanatic  angle  (9Q), 
the  spatial  coherence  length  of  the  atmosphere,  (rQ,  pq)  and 
vertical  profiles  of  the  refractive  turbulence  structure 
parameter,  Cn  .  This  dissertation  deals  with  the  theory  and 
measurement  of  the  isoplanatic  angle  and  the  refractive 
turbulence  structure  parameter.  Currently,  reliable  systems 
exist  for  the  measurement  of  the  spatial  coherence  length  of 
the  atmosphere  [Refs.  1,  2]. 

The  isoplanatic  angle  is  an  angular  measure  of  spatial 
coherence  in  the  atmosphere.  Walters  of  [Ref.  1]  has 
developed  three  generations  of  instruments  (known  as 
i soplanome te r s )  to  measure  the  isoplanatic  angle.  This  work 


has  to  do  in  part  with  the  development  and  test  of  a  fourth 

generation  isoplanometer  having  an  improved  weighting 

function.  This  instrument  is  optimized  for  night  use  and 

provides  temporally  high  resolution  data  (1  second  samples) , 

able  for  the  first  time  to  measure  trends  with  periods  less 

than  five  minutes.  Two  fourth  generation  isoplanometers  are 

currently  operating  at  different  locations. 

Of  the  three  parameters  noted,  the  most  important  is  the 

2 

refractive  turbulence  profile.  Both  the  isoplanatic 

angle,  0Q,  and  spatial  coherence  length,  rQ  are  functions  of 

2  2 
Cn  .  Hence,  knowledge  of  the  Cn  profile  not  only  gives  a 

picture  of  turbulence  with  altitude,  but  also  gives  the 

measures  of  spatial  coherence,  rQ  and  6q.  Unfortunately, 

2 

high  vertical  resolution  profiles  of  Cn  are  difficult  to 

measure.  Several  measurement  techniques  exist,  each  having 

its  strengths  and  weaknesses.  These  techniques  will  be 

discussed  later.  This  research  attempts  remote  measurement 
2 

of  C  by  direct  inversion  of  the  stellar  scintillation 
n 

amplitude  covariance  function.  Mathematically,  the  problem 
is  ill-posed  meaning  that  noise  in  the  data  makes  inversion 
impossible.  Different  regularization  techniques  exist  to 
make  the  problem  tractable.  This  work  is  an  extension  of 
[Refs.  3,  4]  using  the  Tikhonov  regularization  technique 
[Ref.  5]  to  invert  the  covariance  function  for  two  finite 
but  arbitrary  equal  apertures. 


II 


BACKGROUND 


A.  GENERAL  DESCRIPTION  OF  TURBULENCE 

Turbulence  has  certain  characteristics  that  differen¬ 
tiate  it  from  other  flow  patterns.  In  [Ref.  8],  Lumley 
briefly  describes  these  characteristics  as  I  will  do  here. 
The  first  quality  of  turbulence  almost  goes  without  saying; 
turbulence  is  irregular.  This  "quality"  has  a  profound 
effect  on  how  one  approaches  a  mathematical  model  or 
description  of  turbulence.  Fluid  dynamics  for  years  relied 
on  the  application  of  Newton's  laws  to  fluids.  With  the 
Navier-Stokes  equations  in  hand,  one  in  principle  should 
be  able  to  solve  any  problem.  However,  in  the  case  of 
irregular  flows,  the  randomness  makes  this  approach 
impotent  and  we  are  left  with  having  to  use  statistical 
methods . 

Diffusivity  is  a  property  responsible  for  rapidly 
mixing  and  spreading  basic  physical  quantities  like 
momentum,  mass  and  heat.  This  property  is  present  in  all 
turbulent  flows. 


^  rRTor’T*  ~  ~  'a  .v : 


where  v  is  the  velocity,  L  is  a  characteristic  dimension  of 
the  entire  flow  and  v  is  the  viscosity.  As  the  Reynolds 
number  moves  beyond  a  certain  critical  value,  Rec, 

(nominally  Rec  ss  100,  for  pipe  flow)  [Ref.  9],  the  flow 
becomes  unstable  giving  way  to  turbulence.  The  value  of  the 
critical  Reynolds  number  depends  on  the  type  of  flow.  Rec 
for  the  atmosphere  is  approximately  2000  [Ref.  14], 

Another  important  property  of  turbulence  is  that  it  is 
characterized  by  three  dimensional  vorticity  fluctuations. 
The  vorticity  fluctuations  find  their  origin  in  the  velocity 
field.  In  order  to  support  the  random  vorticity  fluctua¬ 
tions,  the  random  velocity  field  must  also  be  three 
dimensional.  Lumley's  example  of  a  two  dimensional 
atmospheric  cyclone  is  worth  noting.  Large  scale  cyclonic 
behavior  responsible  for  weather  patterns  is  itself  not 
turbulent.  This  does  not  mean,  however,  that  the  cyclone  is 
not  made  up  of  smaller  scale  turbulent  bodies. 

Within  every  turbulent  flow  is  a  mechanism  that 
dissipates  energy.  The  dissipation  occurs  at  small  scale 
sizes  and  is  a  function  of  the  viscosity  of  the  fluid.  The 
dissipated  energy  is  internalized  into  the  fluid  typically 
as  heat.  Clearly,  if  one  increases  the  rate  at  which  energy 
is  added  to  a  fluid,  the  dissipation  rate,  e,  must  also  rise 
proportionately  to  conserve  energy.  This  occurs  after  the 
energy  increase  has  reached  the  smallest  scale  sizes. 


-  ■*- 
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The  last  two  properties  in  Lumley's  list  are  rather 
obvious  but  still  should  be  noted.  First,  turbulence  lies 
within  the  continuum  model  of  fluid  mechanics.  The  scale 
sizes  in  turbulence  are  typically  much  larger  than  molecular 
sizes.  Thus,  molecular  diffusion  does  not  significantly 
affect  turbulent  flows.  And  lastly,  turbulent  flows  in 
fluids  are  indeed  flows.  In  principle,  as  I  mentioned 
previously,  the  Navier-Stokes  equations  provide  the  general 
equations  for  fluid  flows.  For  turbulence,  one  is  concerned 
with  the  statistics  of  the  Navier-Stokes  equations. 
Unfortunately,  the  solution  to  this  set  of  equations  is  not 
known  since  the  number  of  unknowns  exceeds  the  number  of 
equations.  This  is  known  as  the  "closure"  problem  of  the 
Navier-Stokes  equations  [Ref.  10]. 

As  was  mentioned  earlier,  turbulence  is  an  irregular  or 
random  process.  Finding  exact  scale  sizes,  dissipation 
rates,  etc.,  is  virtually  impossible.  However,  it  is 
possible  to  extract  order  of  magnitude  values  by  dimensional 
analysis.  This  technique  is  used  in  the  turbulence  field  to 
mathematically  model  the  phenomenon.  Other  techniques 
include  asymptotic  invariance  (the  model  should  behave 
properly  as  the  Reynolds  number  approaches  infinity)  and 
local  invariance  (depending  on  the  time  and  length  scales  in 
the  turbulence,  one  may  be  able  to  assume  that  the 


■  '  .  •  '  . •  >  '  -  -  .*  _»  v  **  >  v  rv*  -  \r:\r  <  «p  r.-  n  *rjr«  rj  »  j 


in  varying  degrees,  find  their  way  into  the  treatment  of 
turbulence . 

All  mechanical  processes  must  have  an  energy  source.  In 
turbulence,  the  source  is  principally  shear  flows  within  the 
mean  flow.  This  is  true  for  turbulence  outside  the 
atmospheric  boundary  layer.  The  boundary  layer  is  formed 
by  the  interaction  of  the  flow  and  the  earth's  surface  given 
the  no  slip  condition  [Ref.  8].  Buoyancy  can  also  be  an 
energy  source.  With  an  energy  source  available,  a  laminar 
flow  can  become  turbulent.  The  laminar  portion  of  the  flow 
is  modeled  by  the  linearized  theories  valid  for  small 
perturbations  much  like  the  classical  or  quantum  harmonic 
oscillator.  The  turbulent  portion  of  the  flow  is  mathe¬ 
matically  treated  using  an  asymptotic  theory  valid  for  high 
Reynolds  numbers.  This  inevitably  leads  to  chaos  in  the 
transition  region. 

In  turbulence  there  are  a  number  of  different  length 
scales.  A  description  of  all  the  different  length  scales  is 
a  key  component  of  similarity  theory  [Ref.  10].  However, 
two  important  ones  are  covered  now.  First,  the  outer  scale, 
typically  noted  as  Lq ,  is  the  size  of  the  turbulence  at 
onset  (on  the  order  of  meters).  The  outer  scale,  given  by 
Lq  =  Re  v/v,  is  a  function  of  the  Reynolds  number,  viscosity 
and  velocity  of  the  fluid  as  mentioned  earlier.  In  the 
outer  scale,  viscosity  does  not  play  a  role  in  the 
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dissipation  of  energy  since  inertial  terms  dominate  the 
equations  of  motion.  However,  at  the  inner  scale  or 
microscale  (on  the  order  of  millimeters) ,  viscosity 
internalizes  the  energy  of  the  turbulence  as  heat.  Clearly, 
from  what  has  already  been  said,  the  key  parameters  at  the 
inner  scale  are  the  dissipation  rate  and  viscosity.  Using 
dimensional  analysis,  Kolmogorov  formed  length,  time  and 
velocity  scales  valid  at  the  microscale  (Ref.  11].  These 
are  called  the  Kolmogorov  microscales  and  their  exact 
expressions  can  be  noted  in  [Ref.  8]. 

Now,  we  can  construct  a  physical  picture  of  the 
formation  and  dissipation  of  turbulence.  The  onset  of 
turbulence  occurs  at  the  outer  scale.  If  the  Reynolds 
number  at  the  outer  scale  is  again  greater  than  the  critical 
Reynolds  number  for  turbulence  formation,  the  turbulence  at 
the  outer  scale  will  again  break  down.  The  formation  of 
continually  smaller  eddies  occurs  until  Re  <  Rec  and  the 
turbulence  is  at  the  inner  scale.  Hence,  energy  from  the 
ordered  flow  moves  through  the  turbulence  "cascade"  to  the 
inner  scale  where  it  is  dissipated  in  random  processes.  The 
region  of  scale  sizes  between  the  inner  and  outer  scale  is 
known  as  the  inertial  subrange.  The  energy  transfer  from 
one  scale  size  to  the  next  is  adiabatic  in  the  inertial 
subrange.  Hence,  the  entropy  increase  occurs  at  the  inner 
scale.  This  view  of  turbulence  is  appealing  since  it  is  so 
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closely  related  to  the  second  law  of  thermodynamics.  We 
might  expect  that  the  inertial  subrange  has  a  well  defined 


spectrum.  It  does,  and  this  will  be  covered  later. 

B.  STATISTICAL  DESCRIPTION  OF  TURBULENCE 

1.  Random  Variables 

Since  turbulence  is  irregular  and  random,  it  is  most 
effectively  described  statistically.  Hence,  variables  like 
the  velocity  v'(t)  ,  become  random  variables  having  a  mean  and 
higher  statistical  moments.  In  many  stochastic  processes, 
the  statistical  moments  do  not  change  with  time  or  are 
stationary.  However,  with  atmospheric  turbulence  this  is 
not  the  case.  The  nonstationarity  of  turbulence  makes  the 
meaning  of  long  term  averages,  etc.,  questionable.  In 
addition,  the  replacement  of  ensemble  averages  by  time 
averages  (ergodic  assumption)  presents  additional 
difficulties.  This  complicates  the  description  of 
turbulence  considerably. 

2.  Homogeneity  and  Isotropy 

Along  with  stationarity,  most  well-behaved  random 
variables  are  homogeneous  and  isotropic.  Homogeneity  means 
that  statistical  quantities  do  not  change  with  a  Galilean 
transformation  of  coordinates.  This  implies  that 
mathematically,  functions  describing  statistics  at  r^ . and  T2 
depend  only  on  the  difference  ?12  =  -  T2.  Isotropy 

implies  a  symmetry  in  rotations  of  r\  and  T-.  Hence,  if  the 


statistics  are  homogeneous  and  isotropic,  functions 
depending  on  and  Y2  can  be  represented  by  17^  -  T2  I . 

Unfortunately,  the  atmosphere  provides  neither 
homogeneous  nor  isotropic  random  variables.  However, 
Tatarski  in  [Ref.  6]  adopts  the  idea  of  local  homogeneity 
and  isotropy.  Over  some  region  R,  comparable  to  the  outer 
scale  Lq,  the  statistics  of  the  random  variables  do  not 
change  with  translations  and  rotations. 

In  a  similar  manner,  Tatarski  handles  the  non¬ 
stationary  problem.  He  notes  that  the  functions  can  be 
stationary  over  time  increments  or  stationary  increments. 
Hence,  it  is  possible  to  tolerate  slow  drifts  in  the  mean 
and  other  moments. 

3.  Structure  Functions 


Another  way  to  cope  with  the  non-stationary  problem 
is  to  define  a  function  in  terms  of  a  difference.  Tatarski 
defines  the  structure  function  for  temperature  as 


dt(T1,t2)  =  <  [T(F2)  -  T(F1)]2  > 


<  >  denotes  an  ensemble  average.  We  should  note  that  the 
structure  function  depends  upon  the  vectors  r^  and  r^.  If 
we  assume  homogeneity  and  isotropy,  this  reduces  to 


Dm(r)  =  <  [T  ( r„)  -  T(r,>r  >  , 


(2.2 


where  r  =  |r*2  -  7^  | .  Also,  in  differencing  T(r^)  and  T(r  j)  / 
effects  of  a  changing  mean  are  removed. 

Kolmogorov  showed  [Refs.  6,  8]  by  dimensional 
analysis  that  the  structure  function  for  temperature 
fluctuations  in  the  inertial  subrange  is  given  by 


DT(r)  =  CT2r2/3  ,  (2.3) 

2 

where  CT  is  the  temperature  structure  parameter.  In  EM 
propagation,  one  is  more  concerned  with  index  of  refraction 
than  temperature.  So,  we  can  define  an  index  of  refraction 
structure  function,  Dn(r),  similar  to  (2.2)  and  (2.3), 

Dn(r)  =  cn2r2/3  ,  (2.4) 

2  . 

where  C  is  the  refractive  turbulence  structure  parameter, 
n 

In  order  to  model  the  turbulence  correctly,  one  must  use 
quantities  that  are  not  affected  by  position  in  the  fluid. 

We  refer  to  such  quantities  as  "passive  additives". 
Temperature  is  not  a  passive  additive  since  it  changes  with 
vertical  displacements.  However,  potential  temperature,  0, 
given  by  e  =  T  -  rz  where  r  is  the  adiabatic  lapse  rate  and 
z  is  altitude,  is  a  passive  additive  since  it  is  corrected 
for  changes  in  altitude. 

Figure  2.1  taken  from  [Ref.  14]  shows  the  behavior 
of  D_(r)  as  a  function  of  log(r).  Clearly,  if  r  <  1  the 


inner  scale,  DT(r)  =  0  and  for  r  >  Lq,  DT(r)  approaches  2 oT 
asymptotically  as  the  fluctuations  become  uncorrelated. 
Hence,  the  structure  function  has  a  well  defined  slope  in 
the  inertial  subrange. 

4 .  Covariance  Functions  and  Spectra 

Along  with  structure  functions,  covariance  functions 
(or  correlation  functions  if  normalized)  and  spectra  are 
useful  in  describing  random  processes.  Conveniently, 
covariance  functions  and  spectra  form  a  Fourier  transform 
pair.  The  covariance  between  two  random  variables  S  and  T 
is  given  by 


Bst  =  <  [T  (rx)  -  <T(r1)>]  [S(r2)  -  <S(r2)>]>,  (2.5a) 

where  the  average  value  of  T,  <T(r)>  has  different  defini¬ 
tions  if  T  is  continuous  or  discrete.  Often  one  uses  the 
autocovariance  function 


Btt  =  <[T(r:)  -  <T(rx)>]  [T(r2)  -  <T(r2)>]>.  (2.5b) 


If  the  random  variable  T  is  homogeneous  or  locally 
homogeneous  B,pT  becomes 


B 


<T(r,  )T(r,)>  -  <T  ( r ,  )  >2  . 


mrp 


(2.6) 


If  we  choose  <T>=0  this  reduces  to 


Btt  =  <T(r1)T(r2)> 


(2.7) 


There  is  a  relationship  between  structure  functions  and 
covariance  functions  given  by,  [Ref.  6], 


DT(r)  =  2[Btt(0)  -  Btt  ( r )  )  .  (2.8) 

As  mentioned  earlier,  the  covariance  function  and 
the  power  spectral  density  form  a  Fourier  transform  pair. 

The  power  spectral  density  gives  one  insight  into  how 
different  frequency  components  contribute  to  the  variance. 
For  one  dimensional  spectra,  the  transform  pair  is 


W(K) 


CO 


(2.9) 


and 


B  ( r ) 


_1 

2ir 


00 


eirK  W(K)  dK 


(2. 10) 


But,  before  these  relations  hold,  we  must  realize  again 
that  the  functions  involved  are  random.  A  random  function 


f(r)  can  be  decomposed  using  Four ier-Stielt jes  integrals 
given  by 


f  (r) 


d*  (K) 


(2.11) 


where  d$(K)  is  a  random  complex  amplitude.  The  Fourier- 
Stieltjes  integral  must  be  used  since  the  random 
fluctuations  are  uncorrelated  over  non-overlapping  intervals 
[Ref.  12].  Hence,  f(r)  does  not  have  a  derivative  and  is 
non-Riemann  integrable. 

With  the  substitution  of  the  Fourier-Stielt jes 
integral,  the  covariance  function  is 


B ( rx+r , rx) 


(K[r1+r] 


K '  r  ] 

1  <d*(K)d**(K')  >  .  (2.12) 


If  the  medium  is  homogeneous,  B(r1  +  r,  r.^)  =  B(r),  since  an 
average  must  not  depend  on  the  location  within  the  medium. 
For  the  double  integral  to  only  depend  on  r,  the  following 
must  be  true. 


We  note 


<d$  (K )  d  <*.  *  (K  '  )  >  =  6  (K-K'  )W(K)  dKdK  '  .  (2.13) 

★ 

that  <d<j>  (K)d<t>  (K')>  is  in  the  form  of  a  mutual 


coherence  function  (MCF) ,  which  represents  the  overlap  of 


* 

<3 ♦  ( K )  and  d*  (K)  .  Since  the  d$'s  are  uncorrelated  except  at 
K  *  K',  we  get  the  Dirac  delta  function  with  the  spectrum 

i 

W{K)  giving  the  amplitude.  Substituting  (2.13)  into  (2.12) 
and  integrating  over  K'  gives  (2.10). 

Since  we  are  dealing  with  fluctuations  of  the  random 
variable  T,  the  mean  value  is  assumed  zero  giving 


2 

CT 


W(K)dK 


(2.14) 


From  (2.8)  and  (2.10)  we  see  immediately  that 


D(r) 


cos(rK)]  W(K)dK 


(2.15) 


since  B(r)  is  an  even  function. 

Tatarski  [Ref.  6]  develops  an  expression  for  W(K) 
based  upon  the  form  of  the  structure  function  for 
turbulence.  The  one-dimensional  spectra  W(K)  is  given  by 

W (K )  =  a  e2/3  K_5/3  ,  (2.16) 


where  a  is  a  constant  and  e  is  the  dissipation  rate. 


Vector  fields  like  the  velocity,  v(r),  are  expressed 


in  three  dimensions.  Hence,  their  spectra  are  three 
dimensional .  We  will  let  *(K)  denote  this  power  spectral 
density.  *  (IT)  is  given  by 


•  (K)  = 


3—  -iK< 
re 


rB(T) 


(2.17) 


The  covariance  function  is  defined  similarly.  B(r)  is  given 
by 


B  ( r )  = 


<2«) 


/// 


d3KelK’r  4 (K) 


(2.18) 


Equation  (2.18)  is  also  developed  from  a  Fourier-Stielt jes 
approach  requiring  homogeneity.  For  spherical  coordinates, 

3—  2 

d  r  =  r  sinedrdedifr  , 


and  the  9  and  4  integrals  in  (2.17)  can  be  completed  giving 


*(K) 


2  it  3K 


B ( r)  sin ( rK) rdr 


(2.19) 


A  similar  relation  for  the  covariance  function  can  also  be 
found.  Tatarski  uses  the  relation 


(2.20 


*(K)  = 


-1  dW (K) 


2irK  dK 


to  evaluate  the  three  dimensional  spectrum  in  terms  of  the 
one  dimensional  spectrum.  This  gives 


»(K)  =  aCn2K_11/3 


(2.  21 


for  the  passive  additive,  index  of  refraction,  where  a  is  a 
constant.  This  is  known  as  the  Kolmogorov  spectrum  (since 
Kolmogorov  did  the  original  work),  valid  in  the  inertial 
subrange.  Spectra  have  also  been  developed  for  the  pro¬ 
duction  and  dissipation  regions  of  turbulence  [Ref.  13]. 


C.  EM  PROPAGATION  THEORY  IN  TURBULENCE 

Clifford  [Ref.  15]  covers  Tatarski's  work  [Ref.  6]  for 
weak  turbulence.  This  treatment  is  analogous  to  any 
application  of  first  order  perturbation  theory.  The 
following  is  Clifford's  treatment  in  a  more  condensed  form. 
1 .  Solving  the  Wave  Equation 

Consider  a  sinusoidal  wave  of  the  form  e~lu,t  where 
the  atmosphere  has  zero  conductivity  and  unit  magnetic 
permeability.  For  these  conditions.  Maxwell's  equation  in 
Gaussian  units  become: 


V  X  H 


-1  kn  E 


t 


(2.24) 


V  •  (n2E)  =  0  ,  (2.25) 

where  k  is  the  wavenumber  of  the  electromagnetic  wave  and  n 
is  the  index  of  refraction.  We  generate  the  wave  equation 
for  If  in  the  usual  way  by  substituting  (2.24)  into  the  curl 
of  (2.23)  giving 

-V2E*  +  V(V  •  E)  =  k2n2F  .  (2.26) 

We  can  expand  (2.25)  and  solve  for  V  •  E.  Substituting  into 
(2.26)  gives 

v2If  +  k2n2If  +  2v(If  •  v  log  n)  =  0  .  (2.27) 

If  we  assume  that  depolarization  effects  are  negligible 
(i.e.  assume  isotropy),  the  last  term  of  (2.27)  can  be 
neglected.  Clearly,  one  need  only  consider  a  single  scalar 
equation,  since  all  the  components  follow  in  the 
mathematics.  Thus,  the  scalar  wave  equation  becomes 

72E  +  k  2n2E  =  0  .  (2.28) 


Clifford  now  applies  the  method  of  small  pertur¬ 


bations.  This  means  that  E  can  be  expressed  as  a  series  of 
decreasing  values. 


E  =  Eg  +  E^  +  E2  +  •  •  * 


(2.29) 


Eg  corresponds  to  the  unscattered  wave,  E^  to  single 

scattering,  E2  to  double  scattering  and  so  on.  In  some 
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situations  the  series  may  diverge,  particularly  if  is 

large  or  the  path  length  long.  This  leads  to  the  phenomenon 
of  saturation.  Clifford  uses  the  first  order  approximation 
giving 


v2E0  +  k2E0  -  0 

and 

v2E1  +  k2E1  +  2k2n1E0  =  0 


(2. 30) 


All  terms  of  order  two  have  been  neglected.  Clifford  like 
Tatarski  assumes  a  unit  amplitude  wave  propagating  in  the  z 
direction.  This  allows  us  to  use  the  standard  Green's 
function  solution  for  E^ (7)  , 


Ei  (r) 


*  zr  (a3r' 


ik|T-T'  i 


[  2k  2n 


<r')eikz'] 


(2.31) 
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For  laser  propagation,  most  of  the  monochromatic  light  will 
be  scattered  into  a  small  forward  cone.  Hence,  backscatter 
is  negligible.  Clifford  now  applies  the  Fresnel 
approximation  to  E^(T)  giving 


E1(r) 


.2  ikz 
k  e 


2  n 


I  d3r '  exp  <  - 
Vol  ( 


2 

/— •  — \  *•  1 


k  [  (T-T' ) 
2  ( z-z  ' ) 


n1  (r ' ) 


(z-z ' ) 


(2. 32) 


where  p*  are  the  coordinates  transverse  to  the  propagation 
direction  z. 

2.  The  Rytov  Approximation 

At  this  point  both  Clifford  and  Tatarski  apply  the 
Rytov  approximation  which  assumes  that  the  solution  of  the 
stochastic  wave  equation  is  of  the  form 


E  =  Aes 


(2.33) 


where  A  is  an  amplitude  and  s  is  a  complex  phase.  The  Rytov 
solution  involves  the  perturbed  wave  equation  (2.30)  divided 
by  Eq.  Tatarski  gives  this  equation  as 


+  k2n2  ( r)  =  7 2  log  E  +  (7  log  E)2  +  k2n2(r)  .  (2.34) 


This  combination  of  solution  and  wave  equation  is  known  as 
the  method  of  smooth  perturbations.  The  name  comes  from  the 


fact  that  we  require  the  perturbations  to  be  small  over  the 
distance  of  a  wavelength.  This  condition  is  less  restrict¬ 
ing  than  those  for  the  method  of  small  perturbations. 

From  (2.33),  the  ratio  of  E/Eq  where  E  =  Eq  +  is 

£ 

i-  -  1  +  =  j-  exp [ i ( S  -  S0)]  .  (2.35) 

E0  E0  A0  ° 

Taking  the  natural  logarithm  of  (2.35)  and  expanding  in  a 
power  series  gives  Clifford's  result. 


*1 

pr-  *  —  +  i  (S  -  Sn)  .  (2.36) 

L0  A0  U 


Thus,  the  Rytov  approximation  has  allowed  us  to  obtain 
amplitude  and  phase  information  from  the  Green's  function 
solution  to  the  wave  equation.  Tatarski  and  others  have 
developed  the  convention 


and 


x  =  tn (A/Aq) 


S 


(2. 37) 


This  is  convenient  since  an  experiment  measures  the  log 
amplitude  x.  It  is  trivial  to  obtain  the  forms  for  A^/Aq 
and  (S  -  Sq)  from  (2.32).  But  it  is  not  so  trivial  to  move 
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on  from  there.  The  problem  is  completing  the  integral  over 
the  random  function  n^(T).  If  one  considers  the  plane  wave 
propagating  in  a  layered  media,  where  at  a  given  z,  there 
exists  statistical  homogeneity,  we  can  expand  n^(r)  in  a 
two-dimensional  Fourier-Stielt jes  integral.  We  also  assume 
that  each  layer  is  uncorrelated  so  that  an  integral  along 
the  path  is  a  sum  of  incoherent  additions.  In  the  notation 
of  Clifford,  nj ( r^  ' )  becomes 


n:  (T  ) 


•/ 


dv  (K ,  z  '  )  e 


iK*  p 


(2. 38) 


where  dv(if,  z')  is  the  complex  Fourier-Stielt  jes  amplitude. 
Substituting  this  into  the  relations  for  A-^/Ag  and  (S  -  SQ) 
gives 


k2  f  f  dz- 

dv (K,z  ' ) 

cos  fk  (7-7* )  *1 

37  J  vz 

(z-z'j  J 

tape 

)  sin  |_ 2  (z-z  ’ ) 

(2.39) 


3 _  2  _ 

where  dr'  =  dz'd  p'.  x(r)  refers  to  the  cosine  of  the 
bracketed  quantity  while  S^  (7)  to  the  sine.  One  can  perform 
the  7*  integrations  leaving  the  Four ier-Stielt jes  integra¬ 
tion  and  the  z'  integration.  Now  we  can  construct  the 
two-dimensional  covariance  functions  exactly  as  before. 
Bx(F,z)  is 


Bx (p,z) 


_ ^  _ 

<X  (  P*2  +  p,  Z)  X  (Pj_ ,  z)  > 


(2.40) 
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Tatarski  finds  the  spectrum  from  (2.40)  to  be 


k  r  2  *| 

Fx  (K,z'-z"=0)  =Fx(K,0)  =  2irk2  J  *  (K)  sin2  — ^z)  d; 


(2.41) 


Fs(K,0)  =  2 irk2  J  *  (K)  cos2 


2  K  (L-z) 


dz  , 


where  L  is  the  path  length  and  *(K)  is  the  three-dimensional 
Kolmogorov  spectrum. 

3.  The  Huvqens-Fresnel  Approach 


As  we  have  seen,  Tatarski  and  Clifford  use  the 
differential  equations  approach  to  solving  the  turbulence 
problem.  Since  we  have  linear  operators,  etc.,  the 
identical  problem  can  be  formulated  in  terms  of  integral 
equations.  This  is  Huygens-Fresnel  theory  applied  to 
propagation  in  turbulence.  Lutomirski  and  Yura  [Refs.  16, 
17]  develop  an  extended  Huygens-Fresnel  theory  by  adding  a 
random  phase  term  to  the  Huygens-Fresnel  integral.  In  fact 
it  is  equivalent  to  applying  the  Rytov  approximation  to  the 
integral  equation.  After  applying  the  paraxial  approxima¬ 
tion,  the  extended  integral  equation  becomes 


.  -ik  / 

(r)  =  ■= —  / 


(ik | r-r ' | )  [*<7)1  2- 

- -  E(?)elHnJ  dr'  , 


(2.42 


where  *  is  a  complex  phase  (equivalent  to  S^,  Yura's 
notation),  the  primes  denote  aperture  coordinates  and  the 
non-primes  denote  image  plane  coordinates.  The  Huygens- 
Fresnel  integral  is  developed  in  standard  texts  like 
[Ref.  18].  We  should  note  that  (2.42)  reduces  to  the 
Green's  function  solution  of  Tatarski  and  Clifford  in  the 
geometrical  optics  limit  (i.e.»  +  ■  k  J* n^(z)dz).  This  is 
not  surprising  since  the  kernal  of  the  Huygens-Fresnel 
integral  is  a  Green's  function. 

Yura  now  finds  the  average  intensity  <I(r)>, 


<1  (r)  >  =  < 
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(2.43) 


The  only  functions  in  (2.43)  that  depend  on  time  are 
the  complex  phases  * '  and  4>".  Hence,  evaluating  the 
ensemble  average,  if  we  assume  ergodicity,  becomes  an  effort 
in  evaluating  the  average  of  the  ^'s.  We  now  change 
variables  to  those  of  Yura  [Ref.  17]  and  Walters  [Ref.  14], 

p+  =  (7'+7")/2  ,  7*  =  7' -7"  ,  7  =  7  . 

With  this  change  of  variables,  the  Fresnel  approximation 


becomes 


(jr-r'l  -  \T  -  T-|)  =  J  (T+  *  T  +  7  *  T~)  ,  (2.44) 


where  z'  =  z"  *  0.  This  gives  a  very  physical  result  for 
the  average  intensity 


2  j 

<I(r)>  =  {-£%)  J  dp7"  ML(T  )MS  (T.z)e  2  ,  (2.45) 


where 


Ml(7  ) 


-  J"e(T*  +  -§— )  E*  CF+  - 


,  i  k  _+  — — . 

?-  <—  »  •“  'u 

— )  e  dp 


(2.46) 


N  e(T,z)  =  <9l*'<T'X)+*m<Tm.Y)l>  # 


(2.47) 


Since  z  >>  r'  or  r",  the  exponential  in  (2.46)  can  be 
approximated  as  unity.  This  makes  ML  the  autocovariance  of 
the  aperture  (or  if  normalized,  the  aperture  mutual 
coherence  function  (MCF) ) .  Mg  is  the  atmospheric  mutual 
coherence  function.  This  allows  us  to  treat  the  propagation 
problem  in  two  distinct  parts;  the  effects  of  the  atmosphere 
and  the  effects  of  the  receiver  optics. 

The  trick  now  is  to  find  the  atmospheric  mutual 
coherence  function.  Yura  assumes  that  the  amplitude 
fluctuation,  x,  and  phase  fluctuation,  s,  that  compose  <l<  are 
Gaussian  variables,  which  implies  that  exp  ( 4» '  +'l* " )  is  log- 
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normally  distributed.  Using  a  property  of  log-normally 
distributed  random  variables  and  some  results  of  Fried's 
work  in  [Ref.  19],  Yura  represents  the  atmospheric  MCF  in 
terms  of  the  wave  structure  function  D(p),  giving 


<e(*'+*"*)>  =  exp(-D(p)/2)  , 


(2.48) 


where  D(p)  =  Dv(p)  +  D_(p).  Now  we  can  apply  the  results  of 
X  s 

Tatarski's  work,  since  we  have  forms  for  F  (K,0)  and  F  (K,0) 

X  5 

which  can  be  related  back  to  their  structure  functions. 

From  this,  D(p)  becomes 


L> 

D(P)  =  3. 44k 2l0~1/3p 2  J  Cn2(z)dz  ;  p  <  P0  (2.49) 


D(p)  =  2.91k2p5/3  j  Cn2(z)  dz  ;  p  >  x/xT  ,  (2.50) 

*^0 


where  Iq  is  the  inner  scale  and  x/x"L  is  the  Fresnel  zone 
size.  Based  upon  (2.49)  and  (2.50),  Yura  represents  (2.48) 


<e  (4-  '+*"*)  >  =  MCF(p)  =  exp  (ZZ)  5/3  ^  (2.51) 


where  pq  is  the  transverse  distance  for  the  MCP  to  decline 
by  e  Yura  has  developed  forms  for  Pq  including  inner  and 
outer  scale  effects  [Ref.  17],  We  should  remember  that  Pq 
is  a  distance  in  the  aperture  of  the  optical  system.  It  is 
nominally  on  the  order  of  a  few  centimeters  during  the  day. 

4 .  The  Modulation  Transfer  Function  Approach 

Fried  in  [Refs.  19,  20]  uses  the  modulation  transfer 
function  (MTF)  rather  than  the  MCF.  This  application  of  the 
MTF  and  MCF  are  closely  related  since  both  represent  the 
effects  of  the  atmosphere  in  different  planes  (i.e., 
aperture  and  image  planes) .  The  MTF  is  the  modulus  of  the 
complex  optical  transfer  function  from  the  linear  systems 
approach  to  optics.  This  points  to  the  fundamental 
difference  between  Fried  and  Yura's  approach.  Fried's  work 
is  done  in  the  coordinates  of  the  image  plane,  i.e.,  spatial 
frequencies,  while  Yura  works  in  the  aperture.  Both  are 
equivalent  by  the  Wiener-Khintchine  theorem  since  a  lens 
produces  the  Fourier  transform  of  the  incident  electric 
field  in  the  image  plane.  The  transformation  between  the 
two  planes  is 

p  ■*  X  R  K  , 

where  R  is  the  focal  length  of  the  optics  and  K  is  the  spatial 
frequency  of  the  scintillations.  This  allows  us  to  move 
between  the  aperture  and  image  planes.  Thus,  the  MTF(K)  is 


■  ■  '7- '.■V'A  l"-  l*-  .V'A  r* 'A  '•  T 


MTF  (K ) 


exp  [-3. 44 


(XRK/r0) 5/3 


(2.52) 


where  Tq  is  the  1/e  distance  of  the  MTF  and  is  given  by 


-5/3 


,2.905, 

6.88  ’ 


k  sec  ( * ) 


OD 

/ 


dzC  (z) 


(2. 53) 


where  *  is  the  zenith  angle  [Ref.  21], 

Explicitly  in  the  MTF  and  implied  in  the  MCF  is  a  low 

pass  filtering  of  the  scintillation  spatial  frequencies  by 

the  atmosphere  and  optics.  This  linear  systems  view  is 

extremely  powerful.  It  allows  us  to  decompose  the  turbulent 

propagation  problem  into  a  product  of  atmospheric  and 

optical  system  filtering  functions.  We  will  use  this  idea 

and  the  Wiener-Khintchine  theorem  later  to  our  advantage. 

5.  Measures  of  Spatial  Coherence 

As  mentioned  in  the  introduction,  Pq  and  Tq 

represent  two  measures  of  the  atmosphere's  ability  to 

maintain  the  spatial  coherence  of  a  propagating  wave. 

Nominal  values  for  r^  are  approximately  3-10  cm  for  ground 

to  space  measurements.  From  (2.53),  we  see  that  turbulence 

close  to  the  aperture  of  an  optical  system  contributes  most 
2  • 

to  rQ  since  Cn  is  largest  in  the  boundary  layer  [Ref.  2). 

Fried  in  [Ref,  21]  derives  the  theoretical  equation 
for  another  measure  of  spatial  coherence  in  the  atmosphere. 
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Isoplanatism  is  the 


This  is  .e  isoplanatic  angle,  9q. 
dependence  of  the  optical  transfer  function  of  a  system  on 
the  position  of  the  source.  Since,  Yura  and  Fried's 
approaches  break  the  propagation  problem  into  MCF's  and 
MTF's  for  the  atmosphere  and  optics,  the  atmospheric 
turbulence  does  not  affect  the  OTF  of  the  optical  system. 

It  affects  the  OTF  of  the  atmospheric  part,  hence  the  term 
anisoplanatism.  This  means  the  atmospheric  OTF  is  dependent 
on  the  particular  propagation  path  through  the  turbulence. 

The  isoplanatic  angle  is  an  angular  measure  of 
spatial  coherence.  If  we  consider  two  paths  through 
turbulence  in  the  atmosphere,  the  isoplanatic  angle  relates 
the  mutual  coherence  e  1  point  between  the  two  paths.  Fried 
defines  the  function  S(1T,^)  given  by 

S  (T,~6)  =  2.905  k  2  J  dvCn2(  |  r  [1-(J)  ]  |  5/3  +  [9V)]5/3 

path 

-  \\  {r[l-(£)]  |2  +  2re[l-(jr)]vc  +  { ©v>2  |  5/6  -  ^|{r[l-(£)]} 

-  2re[l-(£)]vc  +  { e v} 2 | 5/6 )  ,  (2.54) 

where  L  is  the  path  length,  v  is  a  position  along  the  path 
and  c  is  the  cosine  of  the  angle  between  r  and  9.  In  the 
limit  as  r/9  -*■  ® 


(or  two  rays  crossing  at  infinity), 


lim  S  (rVT )  =  2.  905  k  /  dvC  £(0v)J/J  .  (2.5 

r/e  -  -  J  n 

path 

Fried  defines  0Q  as 

e0-5/3  =  2.905  k2  / dvCn2  v5/3  •  (2.5 

path 

Interestingly,  rQ  is  given  by  the  other  limit  of  the 
function  S  (7,TT)  , 


lim  S  ( r  ,  e ) 

r/e  ♦  0 


2.905  k2  f  dvCn2{r[l  -  (^) ] } 5/3 (2.5 
path 


where 


,2.905, 
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k2  /dvC „2ti 

path 
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and  the  (1  -  V/L)  5/^3  term  is  a  spherical  wave  weighting 
function.  Clearly,  we  can  change  variables  from  v  to  z  by 
correcting  for  zenith  dependencies,  i.e.. 


v  *  z  sec  <t>  , 

where  <p  is  the  zenith  angle.  Now  one  can  immediately  see 
the  relationship  between  rQ  and  0Q.  The  isoplanatic  angle 


measurement  looking  up  is  the  same  as  the  rg  or  MTF 
measurement  looking  down  from  space. 


Walters  has  shown  eQ  relates  back  to  pQ  and  the 
Huygens-Fresnel  picture.  The  relationship  is  given  by 

90  =  p0/L  .  (2.59) 

Using  the  Huygens-Fresnel  picture,  one  can  gain  a  physical 
insight  into  rQ  and  0q  measurements.  Consider  a  source  on 
the  earth  propagating  spherical  waves  into  space.  eQ  is  the 
angle  such  that  the  MCF  of  the  field  at  two  different  points 
on  the  wavefront  is  e  1  its  original  value.  This  picture 
can  be  reversed  where  two  points  in  space  emit  spherical 
waves  that  sum  at  a  point  on  the  earth.  This  gives  the 
corresponding  rQ  measurement.  One  of  the  advantages  to  the 
Huygens-Fresnel  approach  is  this  reciprocity. 

Based  upon  the  isoplanatic  angle  weighting  function, 
turbulence  between  8-15  km  is  weighted  heavily.  Hence, 
turbulence  near  the  tropopause  tends  to  make  the  isoplanatic 
angle  small  (i.e.,  3-5  urad).  Large  isoplanatic  angles  are 
in  the  neighborhood  of  13-20  urad.  The  largest  isoplanatic 
angles  (20  urad)  are  very  small,  consequently  adaptive 
optics  concepts  that  compensate  for  the  atmosphere  have  a 
monumental  task  since  the  correction  is  only  good  for  an 
angle  0Q. 

6.  Intensity  Fluctuations 

Tatarski  in  [Ref.  6]  obtains  the  variance  of  the 


intensity  fluctuations  from 


He 


where  B  (0)  is  the  covariance  function  at  p  *  0. 

X 

evaluates  the  integral  over  K  in  B^O)  for  point  sources  and 

point  receivers  giving  results  for  plane  and  spherical 

waves.  In  reality,  we  do  have  point  sources,  stars. 

However,  real  point  receivers  do  not  exist.  Fried  [Ref. 

2 

22],  obtains  a  for  finite  apertures.  He  finds  that  the 
aperture  averaged  variance  is 
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MTF(p)Bx(p) 


(2.61) 


where  D  is  the  diameter  of  the  optics  and  a  is  a  constant 
depending  on  the  type  of  aperture.  Bx(p)  can  be  represented 
by  the  two-dimensional  Hankel  transform  of  the  spectrum 
because  isotropy  imposes  cyclindrical  symmetry. 

Substituting  this  into  (2.61)  gives 


MTF(p)  [ 2 it 


J  JQ  (Kp)Fx (K,0)KdK] .  (2.62) 
0 


The  p  integration  is  trivial  since  we  have  the  Hankel 
transform  of  an  MTF  or  autocorrelation.  This  transform 


yields  the  intensity  spectrum  for  the  finite  aperture 
involved,  I (K) .  Hence,  the  normalized  variance  is 


2ir 


f  KdKFx(K,0)I (K) 
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(2.63) 


where 
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Again  we  see  the  filtering  of  spatial  frequencies  by 

the  atmosphere  in  Fx(K,0)  and  by  the  optics  in  I (K) . 

Aperture  averaging  is  the  filtering  done  by  the  optical 

system.  Equation  (2.63)  has  limits  for  point  apertures  and 

infinite  apertures.  In  [Ref.  6),  Tatarski  shows  that  for 

5  /6 

point  apertures,  (2.63)  behaves  as  z  ,  while  for  infinite 

2 

apertures  it  behaves  as  z  •  This  implies  that  all  other 
apertures  fall  between  these  two  z  dependencies. 

Over  the  years,  research  has  modified  and  shown  the 
limits  of  Fried's  initial  aperture  averaging  work  [Ref.  22]. 
Yura  [Ref.  24]  suggests  a  form  for  the  covariance  function 
for  arbitrary  z.  This  modified  Fried's  results  slightly. 
However,  a  more  serious  problem  was  found  by  Dunphy  and  Kerr 
[Ref.  25]  and  by  Homstad,  et.  al  [Ref.  26],  if  turbulence  is 
very  strong.  Experimental  work  has  failed  to  show  the 
theoretical  aperture  averaging  effects  in  strong  turbulence. 


The  trend  is  towards  less  aperture  averaging  than  expected 
in  strong  turbulence.  Azar,  et.  al.  [Ref.  27]  and  Clifford, 
et.  al.  [Ref.  28]  explain  that  the  covariance  function  for 
strong  turbulence  is  governed  by  two  transverse  scale  sizes. 
The  short  size  determines  the  initial  fast  drop  in  the 
covariance  function,  while  a  second  scale  characterizes  the 
long  tail.  The  long  tail  of  the  covariance  function  reduces 
the  effects  of  aperture  averaging.  The  aperture  averaging 
controversy  will  arise  again  when  we  look  at  isoplanometer 
data . 

7.  Problems  in  the  Theory 

The  theory  for  wave  propagation  presented  so  far  has 

been  for  weak  turbulence  or  single  scattering.  We  have 

neglected  the  effects  of  multiple  scattering.  When  the  path 

lengths  are  long  or  the  turbulence  is  strong,  the  linear 

theory  breaks  down.  Experiments  show  that  for  intensity 
2 

variances  of  a  «  .3,  the  normalized  variance  saturates  to  1 

[Ref.  28].  Multiple  scattering  destroys  the  spatial 

coherence  of  the  wave.  The  linear  theory  assumes  a  coherent 

wave  incident  on  the  perturber  whose  strength  is  dependent 
2 

on  (or  the  temperature  gradient) .  If  the  spatial 

coherence  is  degraded  such  that  it  is  less  than  the  Fresnel 
zone  size,  vTl  ,  the  wave  is  partially  coherent  across  the 


perturber.  This  implies  that  interference  can  no  longer 
completely  modulate  the  wave.  Hence,  the  variance  of  the 


fluctuation  intensity  tends  to  saturate  or  approach  a 
constant.  Other  theories  have  evolved  to  handle  propagation 
in  strong  turbulence  or  multiple  scattering.  Many  of  the 
theories  find  their  basis  in  techniques  developed  for 
quantum  mechanics.  Stroebehn  in  [Ref.  29]  reviews  iterative 
Green's  function  solutions  that  give  an  exact  answer  to  the 
multiple  scattering  problem.  But  the  solution,  in  practice, 
cannot  be  computed.  He  also  summarizes  the  concept  of 
Feynman  diagrams  originally  developed  for  quantum  electro¬ 
dynamics,  and  other  more  sophisticated  theories  developed 
for  strong  turbulence. 
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III.  REMOTE  MEASUREMENTS  OF  THE  ISOPLANATIC  ANGLE 

A.  TECHNIQUE 

Loos  and  Hogge  [Ref.  23]  extend  Fried's  work  of  [Ref. 

22]  by  completing  the  Hankel  transform  of  the  MTF  for  a  "top 
hat"  aperture  function.  Loos  and  Hogge  [Ref.  23]  and 
Walters  [Ref.  30]  realized  that  by  measuring  the  normalized 
variance  of  the  intensity  fluctuations  of  stars  one  could 
obtain  the  isoplanatic  angle  to  within  a  constant.  Walter's 
technique  is  as  follows.  The  normalized  variance  of 
intensity  fluctuations  is 

2  L  ®  2 

-§-  =  4  (4,r2)  (,033k2)  J  Cn2(z)dz  J  KdK  *  (K)  sin2  )  I  (R)  , 

S  0  0 

(3.1) 

where  K  is  the  spatial  frequency  of  the  amplitude 

.  .  .  .  2 
scintillations,  L  is  the  path  length  over  which  Cn  is  non¬ 
zero,  I (K)  is  the  intensity  spectrum  of  the  receiver  optics 
and  ®(K)  is  the  K  portion  of  the  Kolmogorov  spectrum 

for  turbulence  in  the  inertial  subrange.  The  similarity  to 
the  isoplanatic  angle  given  by 


L 

e0_5/3  =  2.  905  k 2 [sec  $]8/3  J  dzCn2(z)z5/3  ,  (3.2) 
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is  obvious 


The  trick  is  to  make  the  weighting  function 


W(z)  , 


W(z)  »  f  KdK«(K)sin2(^2f)I  (K)  ,  (3.3) 

0 

behave  as  z5//3.  Mathematically  this  means  that 

z5/3  «  AW ( z)  ,  (3.4) 

where  A  is  a  proportionality  constant.  Hence,  the 
isoplanatic  angle  becomes 

2 

e0""5/3  -  A'  (-|_)  ,  (3.5) 

S 

where  A'  is  a  new  proportionality  constant  taking  account 
the  four's  and  Pi's,  etc.  of  (3.1). 

Previously  we  noted  that  the  weighting  function  integral 

n  C  /C 

behaves  as  z  for  infinite  apertures  and  z  '  for  point 
apertures  [Ref.  6].  Luckily,  the  isoplanatic  angle 
weighting  function,  z3^3,  lies  within  these  two  limits. 
Hence,  theoretically  it  should  be  possible  to  mimic  this 
behavior . 

The  power  of  this  isoplanatic  angle  measurement  tech- 

2 

nique  is  that  no  explicit  knowledge  of  the  C  profile  is 


2 

required.  Cn  in  the  free  atmosphere  is  difficult  to 
measure  remotely  with  high  vertical  resolution  (as  we  shall 
see  later).  Hence,  integrated  parameters  like  Tq  and  9q 

2  j 

depending  on  Cn  are  usually  measured  independent  of  a  Cn 
profile.  Another  important  advantage  that  naturally  evolves 
from  this  technique  is  the  filtering  done  by  the  aperture 
averaging.  The  averaging  aids  in  suppressing  the  high 
spatial  frequency  scintillations  that  are  most  affected  by 
passage  through  the  atmosphere  and  are  most  liable  to  cause 
saturation. 

B.  THE  FIRST  ISOPLANOMETERS 

Walters  developed  the  first  isoplanometers  to  operate 
during  the  day  and  night  [Ref.  30],  The  optical  receiver  is 
an  Celestron  14  Schmidt-Cassegrain  telescope  obscured  so 
that  a  11  cm  diameter  clear  aperture  remains.  The  theory  at 
the  time  of  instrument  development  was  based  on  a  "top  hat" 
aperture  function.  Walters  found  that  the  11  cm  aperture 
gave  a  good  approximation  to  the  z5^3  weighting  function. 

Since  these  instruments  are  designed  to  operate  at  both 
day  and  night,  they  have  a  very  small  field  of  view  to 
suppress  the  Rayleigh  scattered  background  encountered  in 
daytime  viewing.  Hence,  high  quality  telescope  drives  and 
mounts  are  required  for  accurate  tracking.  In  short,  the 
first  three  generations  of  isoplanometers  require  full-time 
support.  In  [Ref.  30],  Walters  shows  that  these  instruments 
operate  properly  within  the  weak  turbulence  limits. 


One  can  compute  numerically  the  integral  over  K  in 


(3.3)  giving  the  theoretical  weighting  function  for  a 
certain  I  (K) .  In  reality,  the  limits  to  the  K  integral  are 
not  zero  and  infinity.  The  spatial  frequencies  of  the 
scintillations  are  truncated  by  the  inner  and  outer  scales 
of  turbulence.  Hence,  the  upper  and  lower  limits  of 
integration  become 


and 


V 

max 


2  IT 
*0 


K, 


-  2jt 


min 


/ 


where  is  the  inner  scale  and  LQ  is  the  outer  scale.  For 
the  following  computer  simulations,  =  5  mm  and  Lq  =  10  m. 
The  integration  over  K  is  done  by  a  brute  force  application 
of  Simpson's  rule  over  512  points  for  a  given  value  of  z. 

The  z  stepsize  is  log  (Az)  =  .05. 

Figure  3.1  shows  the  error  from  z'*'^  for  a  10  cm  top 
hat  aperture  as  a  function  of  z  for  a  500  nm  electromagnetic 
wave.  I (K )  for  the  10  cm  top  hat  is  the  Airy  pattern, 


where  D  is  the  diameter  of  the  receiver  optics.  The  curve 

in  Figure  3.1  is  normalized  so  that  the  integrated  error 

between  2  and  20  km  is  zero.  This  ensures  that  the  error  is 
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not  consistently  above  or  below  z  giving  erroneous 

results  when  this  weighting  function  is  integrated  into  a 

measurement  system.  The  normalization  interval  2  to  20  km 

2 

results  from  the  fact  that  turbulence,  hence  Cn  ,  is  large 

enough  in  that  region  to  contribute  to  the  isoplanatic  angle 

5/3  2 

integral  given  the  z  dependence.  Clearly,  is  largest 

in  the  boundary  layer  [Ref.  2],  however  the  weighting  is 

very  small  there  causing  little  change  in  the  isoplanatic 

2 

angle.  Actual  profiles  of  Cn  ,  like  that  provided  by  Dr.  E. 
Good  from  the  Air  Force  Geophysics  Laboratory,  Hanscom  AFB , 
MA  in  Figure  3.2,  bear  out  this  assumption.  The  particular 
profile  shown  will  be  used  extensively  in  Chapters  IV  and 
VIII  of  this  dissertation.  A  more  complete  description  of 
the  profile  will  be  provided  at  that  time. 

In  [Ref.  24],  Loos  and  Hogge  plot  the  two 
theoretical  limits  for  point  and  infinite  apertures  as  a 
function  of  z.  They  also  plot  the  aperture  averaged 
weighting  function  for  a  35  cm  telescope.  They  note  that 

2  5/3 

the  large  aperture  is  on  the  z  side  of  z  .  In  Figure  3.1 

2 

we  see  two  different  slopes  present.  If  W(z)  =  z  ,  the 
slope  would  be  +1/3,  while  if  W(z)  =  ,  the  slope  would 


be  -5/6.  At  low  z,  the  slope  in  Figure  3.1  is  0.26,  while 


as 


at  high  z  the  slope  is  -0.63.  As  the  diameter  gets  larger 

the  slope  of  the  low  z  portion  moves  closer  to  1/3. 

Clearly,  these  figures  only  show  the  general  influence  of 

the  limits.  Figure  3.3  shows  the  movement  of  the  peak 

towards  higher  z  with  increasing  aperture  size.  Thus,  we 
2 

see  more  z  behavior  as  expected. 

The  question  arises  as  to  why,  with  increasing 

2 

aperture  size,  one  sees  more  z  behavior.  The  answer  is 

straightforward.  From  [Ref.  23],  we  see  that  as  D,  the 

diameter  of  the  optics  goes  to  infinity,  I(K)  =  6(K-0).  As 

finite  D  increases,  the  width  of  the  Airy  disk  narrows  due 

to  the  Hankel  transform.  Eventually,  the  Airy  disk  becomes 

2 

a  good  approximation  to  the  6-function,  hence,  the  z 
behavior . 

Figure  3.4  shows  the  wavelength  dependence  of  the 
weighting  function.  To  first  order,  the  normalized  variance 

given  by  (3.1)  is  wavelength  independent.  This  can  be  shown 

2  2  . 
by  expanding  sin  (K  z/2k)  in  a  Maclaurin  series  giving 


c,  ,  .  2,K2z, 

f(z)  =  sin  (-^-) 


4  2 

- 

4kz 


(3.7) 


2 

The  k  cancel  in  (3.1)  removing  the  wavelength  dependence. 

2 

We  also  see  from  (3.7)  how  the  z  weighting  function 
mentioned  above  comes  about.  The  weighting  function,  (3.3), 


5 


is  itself  not  wavelength  independent,  hence  we  see  the 
behavior  persent  in  Figure  3.4.  As  the  wavelength  increases 
the  weighting  function  moves  laterally  towards  lower  z. 

The  calibration  constant.  A',  from  (3.5),  for  the 
10  cm  top  hat  aperture  is  nominally  1.0  x  10^®.  We  would 
expect  that  this  constant  should  get  larger  as  more  aperture 
averaging  is  included  since  the  normalized  variance  will 
decrease. 

2.  Improved  Weighting  Functions 

The  error  from  z5^3  in  Figure  3.1  at  100  m  is  about 

a  factor  of  three.  In  the  weighting  function  given  by 

(3.3),  the  only  variable  that  depends  on  the  measurement 

system  is  I(K).  The  remainder  is  due  to  the  amplitude 

fluctuation  spectra  Fx(K,0).  For  some  time,  different 

groups  have  been  modifying  telescope  apertures  or  image 

2 

planes  to  measure  Cn  profiles  [Refs.  31,  32,  33].  The 
modification  of  aperture  functions  by  masks,  etc.,  is  called 
aperture  apodization.  Since  the  z5//3  weighting  function 
lies  between  the  theoretical  limits,  we  hope  that  through 
aperture  apodization  we  might  improve  the  isoplanometer 
weighting  function. 

To  obtain  intensity  spectra  from  different  aperture 
functions,  one  can  Hankel  transform  the  MTF  or  autocorrel¬ 
ation  of  the  aperture  or  Hankel  transform  the  aperture 
function  and  square  the  resultant  amplitude  spectra.  The 
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results  are  equivalent  by  the  Wiener-Khintchine  theorem. 
Bateman  [Ref.  34]  compiled  many  different  Hankel  transform 
pairs.  Using  this  list  and  different  combinations  of  Bessel 
functions,  I  began  constructing  different  weighting 
functions.  Appendix  A  is  a  synopsis  of  the  some  apertures, 
spectra  and  resultant  weighting  functions.  We  found  that 
combinations  of  Bessel  functions  are  most  desirable  for 
spectra,  since  they  represent  apertures  that  can  be 
implemented  easily  on  a  telescope. 

Before  moving  on  to  the  aperture  function  that 

approximates  the  z5//^  quite  well,  it  is  interesting  to  look 
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at  the  error  from  z  '  produced  by  a  simple  annulus 
aperture.  Figure  3.5  is  an  error  plot  for  an  annulus  with  a 
10  cm  outer  diameter  and  a  4  cm  obscuration  diameter.  The 
two  peaks  correspond  to  the  influence  of  the  inner  (lower  z 
peak)  and  outer  (higher  z  peak)  diameters  of  the  two  Bessel 
functions  involved.  This  leads  one  to  believe  that  by  using 
combinations  of  Bessel  functions  the  dependence  could 

be  mimicked.  Figure  3.6  is  the  error  for  a  Celestron  C8 
telescope  aperture.  It  is  interesting  to  note  that  the 
inner  diameter  affects  the  weighting  function  at  high  z, 
while  the  outer  diameter  affects  performance  at  low  z.  We 
saw  this  before  with  the  two  theoretical  limits  of  the 
weighting  functions. 

After  some  effort,  an  aperture  function  was 
identified  for  a  8"  Schmidt-Cassegra in  telescope  that 


produces  a  weighting  function  that  approximates  the  z 
dependence.  The  aperture  is  the  double  annulus  shown  in 
Figure  3.7  producing  the  error  from  z5//3  shown  in  Figure  3.8 
for  a  500  nm  electromagnetic  wave.  Over  the  altitude  range 
critical  for  isoplanatic  angle  measurements  (2-20  km) ,  this 
weighting  function  does  very  well. 

The  weighting  function  for  the  double  annulus 
aperture  function  is  given  by 


W(z)  = 
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where  e^,  e2,  are  constants  that  ratio  the  inner 
diameters  to  D,  the  outer  diameter  of  the  optics.  *(K)  is 
the  K_11,/3  dependence  of  the  Kolmogorov  spectrum.  Figure 
3.9  shows  a  comparison  of  the  amplitude  spectra  squared 
(intensity  spectra)  for  the  20  cm  double  annulus  and  a  20  cm 
top  hat  aperture.  It  is  interesting  to  note  that  the 
increased  higher  spatial  frequencies  present  in  the  double 
annulus  spectrum  are  responsible  for  the  improvement  in  the 
weighting  function. 
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3.  Double  Annulus  Calibration 


The  calibration  constant,  A',  for  the  20  cm  double 

annulus  aperture  function  is  approximately  a  factor  of  four 

larger  than  the  10  cm  top  hat.  This  is  due  to  the  increased 

aperture  averaging.  The  double  annulus  weighting  function 

shows  the  same  wavelength  dependence  as  the  top  hat.  Figure 

3.10.  An  increase  in  the  wavelength  of  the  EM  radiation 

shifts  the  weighting  function  to  lower  z.  One  interesting 

point  is  that  the  normalized  variance  has  no  wavelength 

dependence  as  we  have  shown.  However,  the  isoplanatic  angle 

certainly  does  (3.2).  Clearly,  all  the  wavelength 

dependence  for  the  isoplanatic  angle  is  in  A'  and  it  should 

2 

behave,  to  first  order  as  k  .  Figure  3.11  shows  the  A' 

wavelength  dependence.  At  wavelengths  greater  than  800  nm, 

2 

A'  behaves  as  k  but  diverges  below  800  nm.  Clearly,  this 
is  important  depending  on  the  temperature  of  the  stars  used 
for  isoplanatic  angle  measurement. 

Based  upon  the  factor  of  four  increase  in  the 
calibration  constant,  one  might  expect  that  an  instrument 
with  this  weighting  function  would  be  able  to  withstand  a 
factor  of  four  more  turbulent  intensity  before  saturating. 
However,  (Refs.  25,  26,  27]  show  that  the  aperture  averaging 
effects  may  be  degraded  in  strong  turbulence  making  the  10 
cm  and  20  cm  devices  saturate  at  the  same  relative 
normalized  variance.  This  issue  will  be  resolved  later  by 
experiment. 
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Figure  3.10  Wavelength  Dependence  of  the  Error  from  z  for  a  Double  Annulus 


The  20  cm  double  annulus  has  approximately  twice  the 
area  of  the  10  cm  top  hat.  This  implies  that  the  telescope 


will  collect  twice  the  photons  increasing  the  average 
intensity  by  a  factor  of  two.  The  variance  of  the  shot 
noise  current  will  also  increase  by  a  factor  of  two. 
Mathematically,  the  isoplanatic  angle  for  the  10  cm  system 
is 


2  2 
where  a  is  the  variance  of  the  signal  and  <j  is  the 
s  n 

variance  of  the  noise.  For  the  20  cm  double  annulus  system, 
this  same  equation  becomes 


2an 
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The  proportionality  constant,  A’,  for  the  20  cm  system  is 
approximately  a  factor  of  four  larger  than  for  the  10  cm 
system.  In  practice,  we  measure  the  shot  noise  experi¬ 
mentally  and  remove  it  from  the  variance  in  software. 

Hence,  its  effect  is  negligible  at  night.  The  10  cm 
isoplanometer  uses  a  pellicle  to  split  the  light  beam  in  the 
detector  subassembly.  This  reduces  the  average  intensity  by 
a  factor  of  two.  The  20  cm  system,  without  the  pellicle,  is 
equivalent  to  the  10  cm  system. 


4.  Double  Annulus  Improvements 


The  single  greatest  advantage  the  double  annulus 
aperture  function  isoplanometer  would  have  over  the  previous 
systems  is  that  the  35  cm  telescope  is  no  longer  required. 
This  fact  alone  means  saving  approximately  $8000.00  per 
system,  not  to  mention  that  an  isoplanometer  based  on  an  8" 
telescope  can  be  set  up  and  operated  by  a  single  individual. 
The  weighting  function  itself  is  a  significant  improvement 
over  that  of  the  top  hat.  We  will  see  in  the  next  chapter 
just  what  this  improvement  means  in  terms  of  isoplanatic 
angle  measurement. 


IV.  APODIZED  ISOPLANOMETER  SENSITIVITY  ANALYSIS 

Two  questions  might  be  asked  based  on  the  double  annulus 

weighting  function  of  Figure  3.8.  First,  at  100  m  there 

2 

remains  a  factor  of  1.75  error.  C  at  low  altitudes  is 

n 

large,  so  deviations  from  z may  cause  significant  errors 
in  the  isoplanatic  angle  measurement.  Second,  in  the 
numerical  calculations  of  the  weighting  functions,  a  fixed 
inner  scale  was  used  for  all  z.  Clearly,  the  inner  scale  is 
a  function  of  altitude  which  may  drastically  change  the 
weighting  function.  These  two  issues  are  addressed  in  this 
chapter. 

A.  WEIGHTING  FUNCTION  PERFORMANCE  IN  STRONG  TURBULENCE 
Consider  an  isoplanometer  viewing  a  star  directly 

2 

through  a  convective  plume.  Clearly,  the  values  of  Cn  will 
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be  very  large;  on  the  order  of  10  m  near  the  ground. 

From  Figure  3.8,  one  can  immediately  see  that  the  weighting 

function  for  the  double  annulus  has  its  greatest  error  in 
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this  region.  Even  though  the  z  weighting  is  low  near  the 

2 

earth,  the  Cn  values  in  strong  turbulence  may  be  large 
enough  to  cause  significant  error  in  isoplanatic  angle 
measurements.  Based  on  this  concern,  I  performed  a 


sensitivity  analysis  on  both  the  10  cm  top  hat  and  double 
annulus  aperture  weighting  functions. 


The  C  Profile 


_ n _ 

The  most  realistic  way  to  test  the  weighting 

2 

functions  is  to  use  real  Cn  data.  Figure  3.2  shows  Good's 
2 

C  data  between  500  m  and  40  km.  Good  obtained  this 
n 

profile  by  differential  temperature  measurement  between  two 
balloon-borne  microthermal  probes.  From  these  differences 

one  can  construct  the  temperature  structure  function  DT  and 

2  2  2 
obtain  C  as  a  function  of  z.  C  and  C_  are  related  by 
n  n  t 

the  expression,  [Ref*  6], 

Cn2  =  ct2<79p/t2  x  I9"6)2,  (4.1) 

where  P  is  the  pressure  in  millibars  and  T  is  the 

2 

temperature  in  Kelvin.  Direct  measurement  of  Cn  in  this 

manner  yields  very  high  vertical  resolution,  on  the  order  of 

tens  of  meters.  AFGL  averaged  raw  data  to  give  the  500  m 

resolution  in  Figure  3.2. 

For  purposes  of  this  analysis,  the  region  of 

greatest  interest  lies  below  500  m.  Walters  [Ref.  2],  plots 
2 

C  versus  altitude  for  1  to  1000  m.  Walters  measured  this 
n 


particular  data  set  at  midday  in  a  desert  location  and  it 

represents  very  strong  turbulence  conditions.  To  account 

for  convective  plumes  in  low  altitude  turbulence,  I  doubled 
2  2 

the  Walters  C  values.  This  gives  one  C„  profile  from  1  m 
n  3  n 

to  40  km. 
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2.  The  Analysis 

For  convenience,  I  segmented  the  profile  into  four 
regions,  1  to  10  m,  10  to  100  m,  100  to  500  m,  and  500  m  to 
40  km.  Over  each  of  these  regions,  I  numerically  evaluated 
three  integrals, 


Theoretical 


Top  Hat  Aperture 


Double  Annulus  Aperture 


Cn2(z)z5/3dz  , 

Cn2(z)WTH(z)dz  ' 
2(z)WDA(z)dz  . 


(4.2) 


Based  on  Figures  3.1  and  3.8,  we  would  expect  the  weighting 

function  integrals  to  be  smaller  than  theoretical  at  low  z. 

Table  1  shows  the  results  of  this  analysis.  Table  2 
contains  the  relative  errors  for  each  segment  and  the  total 

path.  The  integrated  values  in  Table  1  show  that  even  the 

strongest  turbulence  at  low  z  accounts  for  only  a  small 
percentage  (about  2%)  of  the  integrated  path  effect. 

Clearly,  the  errors  over  the  total  path  for  both  the  top  hat 
aperture  and  the  double  annulus  are  within  the  limits  of  the 
theory.  As  one  might  expect,  this  analysis  is  extremely 
sensitive  to  the  choice  of  A,  a  proportionality  constant. 


where 


TABLE  1 


CONTRIBUTIONS  TO  THE  ISOPLANATIC  ANGLE  INTEGRAL  FOR 
STRONG  TURBULENCE  ANALYSIS 


Theoretical 

Top  1 

Hat 

Double  Annulus 

1-10  m 

2.  61 

X 

10'11 

3.44 

X 

io'12 

5.  80 

X 

CM 

H 

1 

O 

H 

10-100  m 

8.19 

X 

ID'10 

2.28 

X 

io’10 

3.77 

X 

io"10 

100-500  m 

1. 18 

X 

io-8 

5.  50 

X 

IO'8 

8.67 

X 

10~9 

500  m-40  km 

3.89 

X 

10-7 

3.88 

X 

10-7 

3.91 

X 

10~7 

TABLE  2 

RELATIVE 

ERROR  FROM  THEORETICAL 
TURBULENCE  ANALYSIS 

FOR  STRONG 

Top  Hat 

Double  Annulus 

1-10  m 

r- 

00 

• 

CO 

r- 

• 

10-100  m 

.72 

.54 

100-500  m 

.69 

.  26 

500  m-40  km 

1 

O 

X 

r- 

6  x  10-3 

1  m-40  km 

.02 

3  x  10-3 
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In  building  an  isoplanometer ,  the  better  one  can  define  A  or 
A',  the  better  the  instrument  should  operate. 

This  analysis  shows  that  an  isoplanometer  should 
perform  well,  based  on  its  weighting  function  in  strong 
turbulence.  However,  we  must  note  that  in  very  strong 
turbulence,  saturation  has  already  reduced  the  reliability 
of  the  isoplanatic  angle  measurements. 

B.  INNER  SCALE  EFFECTS 

The  inner  scale,  n,  is  given  by 


n 


v3  1/4 
<7“> 


(4.4) 


where  v  is  the  kinematic  viscosity  and  e  is  the  energy 
dissipation  rate.  Clearly,  the  kinematic  viscosity  and  the 
dissipation  rate  are  both  functions  of  altitude.  In  com¬ 
pleting  the  weighting  function  integral  (3.3)  numerically,  I 
used  the  same  inner  scale  for  all  values  of  z.  This  could 
mean  that  the  weighting  function  does  not  accurately  model 
the  z5/3  behavior  required  for  isoplanatic  angle 
measurement . 

1 .  Constructing  the  Inner  Scale  as  a  Function  of 
Altitude 

From  (4.4),  we  see  that  the  inner  scale  depends 


strongly  on  the  kinematic  viscosity  v,  given  by 


where  e  is  the  viscosity  coefficient  and  p  is  mass  density. 

The  U.S.  Standard  Atmosphere  [Ref.  35]  has  vertical  profiles 

of  the  kinematic  viscosity.  Hence,  this  portion  of  the 

inner  scale  is  well  understood. 

The  dissipation  rate  is  another  story  entirely. 

Turbulence  in  the  free  atmosphere  is  intermittent, 

temporally  and  spatially,  unlike  boundary  layer  turbulence. 

This  is  compounded  by  the  fact  that  measurements  in  the  free 

atmosphere  are  difficult  to  make.  Using  [Refs.  36-41],  I 

constructed  what  might  be  called  a  "nominal"  dissipation 

rate  profile.  This  profile,  by  no  means,  represents  the 

output  of  any  theoretical  model  on  my  part.  It  is  a 

synopsis  of  the  data  and  model  results  given  in  [Refs. 

36-41].  The  "nominal"  dissipation  rate  profile  is  shown  in 

Figure  4.1.  Some  controversy  exists  as  to  the  behavior  of 

the  dissipation  rate  at  higher  altitudes  [Ref.  38].  It  is 

generally  agreed  that  e  can  reach  values  comparable  to  that 

observed  in  the  surface  layer.  My  profile  shows  light- 
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moderate  turbulence  near  the  surface,  e  =  40  cm  /sec  (or 

2  3 

log  e  =  1.6)  at  100  m  and  decreases  to  .38  cm  /sec  (or 

log  e  =  -0.42)  at  55  km.  As  a  point  of  reference,  we  might 

note  that  the  dissipation  rates  in  convective  storms  are  on 
3  4  2  3 

the  order  of  10  -10  cm  /sec  [Ref.  40],  The  units  on 


dissipation  rate  indicate  a  "specific"  dissipation  rate 
where  system  mass  is  removed. 

Based  on  the  Standard  Atmosphere  kinematic  viscosity 
profile  and  my  "nominal"  dissipation  rate  profile,  one  can 
now  construct  the  behavior  of  the  inner  scale  with  altitude. 
Figure  4.2  shows  the  inner  scale  profile.  Clearly,  the 
inner  scale  changes  dramatically  with  altitude.  At  55  km, 
n  is  on  the  order  of  a  meter.  As  expected,  the  inner  scale 
follows  the  kinematic  viscosity  since  kinematic  viscosity 
increases  three  orders  of  magnitude  between  100  m  to  50  km. 

2.  Weighting  Function  Dependence  on  the  Inner  Scale 

A  crude  approximation  to  the  inner  scale  effect  can 
be  obtained  by  increasing  the  inner  scale  for  all  z  while 
numerically  evaluating  the  weighting  function  integral. 
Figure  4.3a  shows  how  the  double  annulus  aperture  weighting 
function  changes  with  an  inner  scale  of  10  cm.  The 
deviation  in  W(z)  caused  by  the  change  in  inner  scale  from  5 
mm  to  10  cm  is  not  catastrophic,  however,  a  more  complete 
analysis  needs  to  be  done. 

Since  the  upper  limit  of  the  K  integration  depends 
on  the  inner  scale,  we  see  that  increasing  inner  scale  size 
corresponds  to  lowering  the  upper  limit  of  integration.  We 
expect  that  W(z)  will  not  change  significantly  at  lower  z, 
due  to  the  modest  increase  in  the  inner  scale.  However,  as 
the  inner  scale  increases,  the  upper  limit  of  integration 


DOUBLE  RNNULUS  RPERTURE  FUNC 
ER  DIRM= .0750  (M);  OUTER  DIRM* 


will  decrease.  At  some  critical  Kmax,  the  upper  limit  will 
begin  eliminating  integrand  producing  a  dramatic  change  in 
the  weighting  function. 

Figure  4.3b  shows  the  error  from  z**^  for  the  double 

annulus  aperture  weighting  function  with  the  z  dependence  in 

the  upper  limit  of  integration.  As  expected,  we  see  little 

change  at  lower  altitudes.  At  40  km,  the  inner  scale  is  now 

large  enough  to  begin  removing  part  of  the  integrand.  The 
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deviation  from  z  begins  when  the  inner  scale  reaches  15 
cm,  and  by  25  cm,  the  error  from  z5/^  is  a  factor  of  three. 

The  inner  scale  behavior  of  the  weighting  function 

is  not  a  problem  in  a  measurement  system,  since  values  of 

2  -22  -2/3 

Cn  at  40  km  are  approximately  10  m  7  .  Disturbances 

to  the  electromagnetic  wave  propagating  through  the  medium 

2 

occur  at  lower  z  where  C  is  larger. 

n 

C.  SENSITIVITY  ANALYSIS  CONCLUSIONS 

Based  on  the  results  of  sections  A  and  B,  the 
isoplanatic  angle  weighting  function  appears  very  robust. 

The  only  real  theoretical  problem,  in  terms  of  measurement, 
is  saturation  of  the  normalized  variance.  One  might  ask 
about  outer  scale  effects.  Nominally,  the  outer  scale, 
except  very  near  the  surface,  remains  on  the  order  of  10  m 
[Ref.  42].  This  can  be  seen  from  the  Reynolds  number.  The 
Reynolds  number  is,  (as  shown  in  the  Background  section). 


CTION 
= . 1000 


bh  Al 


where  v  is  the  velocity,  L  is  a  characteristic  size  of  the 
flow  and  v  is  the  kinematic  viscosity.  The  onset  of 
turbulence  occurs  at  Reynolds  numbers  of  2000  or  greater, 
in  the  free  atmosphere  remains  about  the  same  and  v 
increases  with  altitude.  Thus,  velocity  is  the  parameter 
responsible  for  large  Re.  Since  L  does  not  vary 
considerably,  LQ,  the  next  scale  size  smaller,  should  not 
vary  considerably  either.  Near  the  earth,  L  has  an 
immediate  fixed  boundary  condition  that  reduces  its  size. 


Currently,  many  DoD  agencies  require  isoplanatic  angle 
measurements  at  night.  Based  on  this  need,  we  fabricated  a 
fourth  generation  isoplanometer .  What  follows  is  a  brief 
description  of  the  hardware,  detailed  analysis  of  instrument 
qualification  and  a  comparison  with  the  10  cm  "top  hat" 
isoplanometer. 

A.  HARDWARE 

1.  Telescope 

The  double  annulus  aperture  function  of  Figure  3.8, 
was  developed  for  use  on  an  8"  Schmidt-Cassegrain  telescope. 
The  Celestron  C8  is  a  reasonably  priced  and  readily 
available  Schmidt-Cassegrain  telescope  with  a  high  quality 
worm-gear  drive  for  star  tracking.  The  choice  of  Bessel 
functions  as  a  spectral  basis  set  has  the  advantage  of 
simple  obscuration  implementation.  The  aperture  obscuration 
is  two  disks  made  of  flat  black  rubber  adhered  to  the 
Celestron  corrector  plate  with  rubber  cement.  This  allows 
removal  at  a  later  date  without  damage  to  the  optics  or 
coatings.  With  the  exception  of  the  obscuration,  all  of  the 
equipment  associated  with  the  telescope  (tripod,  mount, 
etc.)  is  furnished  by  Celestron. 
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2.  Detector 


The  detector  subassembly  is  a  modification  of  the 
already  existing  10  cm  "top  hat"  isoplanometer  detector 
developed  by  Walters  [Ref.  30].  Walters  engineered  and 
implemented  the  modifications  to  the  fourth  generation 
detector  subsystem,  associated  signal  conditioning  hardware 
and  software.  The  detector  subassembly  screws  into  the  rear 
of  the  telescope.  Incident  radiation  is  directed  over  one 
of  two  paths.  To  acquire  stars  and  focus  the  telescope,  the 
light  enters  a  wide  angle  eyepiece.  Otherwise,  the  light 
passes  through  a  pupil  (on  the  order  of  millimeters),  and  is 
incident  on  a  photomultiplier  tube  that  senses  the  intensity 
fluctuations  (stellar  scintillations) .  The  image  is 
slightly  defocused  to  remove  the  effects  of  any  inhomogen¬ 
eities  on  the  surface  of  the  photomultiplier.  The  photo¬ 
multiplier  is  chosen  to  have  good  red  response,  which  causes 
a  slight  improvement  in  the  weighting  function  at  low  alti¬ 
tudes,  Figure  3.10  and  suppresses  any  Rayleigh  backgound. 

We  should  note  that  the  pupil  in  this  instrument  can 
be  much  larger  than  the  previous  isoplanometer,  since  it  is 
intended  for  night  use  only.  The  10  cm  isoplanometers  use  a 
chopper  to  make  reference  background  measurements  necessary 
Ijr  daytime  operation.  This  feature  is  not  required  on  the 
fourth  generation  system  for  the  same  reason  as  above.  A 
schematic  diagram  of  the  detector  subassembly  is  shown  in 
Figure  5.1. 
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Fourth  Generation  Isoplanometer  Detector  Schematic 


3.  Signal  Conditioning  and  Isoplanatic  Angle  Extraction 

The  first  fourth  generation  isoplanatic  angle 
measurement  system  (A  system) ,  uses  a  Hewlett-Packard 
HP  3421  Data  Acquisition  unit  to  sample  the  output  of  a 
signal  conditioning  unit.  The  signal  conditioning  unit 
performs  two  functions.  It  low  pass  filters  the  incoming 
signal  below  500  Hz  and  calculates  the  mean  and  rms  values 
of  the  input  voltage  data  from  the  photomultiplier.  The 
Data  Acquisition  Unit,  with  an  internal  voltmeter,  detects 
the  voltage  on  each  channel  once/second  and  relays  the 
information  to  a  Hewlett-Packard  217  Computer. 

Once  the  data  is  in  the  computer,  the  normalized 
variance  of  the  signal  is  calculated  from  the  mean  and  rms 
values  of  the  voltage.  Obtaining  the  isoplanatic  angle  is 
then  a  simple  matter  of  applying  the  constant  A'  to  the 
normalized  variance.  The  average  intensity  (in  volts)  and 
isoplanatic  angle  (in  yrads)  are  plotted  on  a  graphics 
monitor.  A  real-time  plot  of  the  average  intensity  can  be 
very  helpful  during  data  acquisition  in  determining  the 
occurrence  of  clouds  and  monitoring  stellar  tracking  errors. 
The  isoplanatic  angle  is  very  robust  to  drifts  in  the 
average  intensity  since  we  normalize  the  variance. 


On  later  systems  (B  systems) ,  the  HP  3421  is 
replaced  by  an  Infotek  A-D  converter  that  resides  within  the 
HP  computer.  The  computer  samples  the  A-D  converter  so  that 


one  isoplanatic  angle  measurement  is  made  per  second.  The 
computer  does  not  store  all  the  isoplanatic  angle  data. 
Instead,  the  computer  calculates  ten  second  averages  and 
stores  the  averages  every  15  minutes. 

The  fourth  generation  isoplanometer  uses  software 
that  has  evolved  and  become  more  "friendly"  with  each  field 
experiment.  Currently,  the  zenith  angle  dependence  of  the 
isoplanometer  angle  is  removed  within  the  main  program. 


B.  ZENITH  ANGLE  QUALIFICATION  TESTS 

It  is  an  interesting  problem  to  know  if  your  instrument 

is  really  measuring  the  isoplanatic  angle.  One  could  fly 

2 

balloons  and  obtain  a  direct  measurement  of  the  Cn  profile 
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and  integrate  with  z  weighting,  however,  Walters  [Ref. 
30] ,  has  developed  a  much  less  costly  and  simpler  technique 
The  isoplanatic  angle  is 


-5/3 


=  [sec  tj)] 


8/3 


CTs2 

(“5") 


(5.1 


including  the  zenith  angle  dependence.  For  increasing  path 

length  (higher  zenith  angles),  the  normalized  variance  will 

increase  since  the  light  is  encountering  more  integrated 

2  2 

turbulence.  Hence,  a  plot  of  os  /S  versus  log  (sec  <j> ) 
should  have  an  8/3  slope  if  the  instrument  is  truly 
measuring  the  isoplanatic  angle  to  within  a  multiplicative 


constant . 


Walters  used  this  technique  to  qualify  the  10  cm  "top 
hat"  isoplanometer  (really  11  cm)  [Ref.  30].  The  following 
is  an  application  of  the  above  technique  to  the  fourth 
generation  system,  but  first  a  caution.  Clearly,  if  the 
atmosphere  is  horizontally  homogeneous  (i.e.,  composed  of 
onionskin  layers) ,  the  technique  should  be  effective. 
However,  if  the  atmosphere  is  irratic  and  the  turbulence  is 
not  uniform  over  the  night  sky,  this  technique  will  not  be 
as  reliable. 

We  (D.  L.  Walters  and  myself)  obtained  qualification 
measurements  for  the  fourth  generation  system  at  two 
different  locations.  The  first  data  set  is  from  McDonald 
Observatory  in  western  Texas.  The  night  was  ideal  for 
instrument  qualification  using  the  zenith  angle  technique. 
The  isoplanatic  angles  were  on  the  order  of  15y  rad.  In 
Figure  5.2  each  point  represents  one  10  sec  average  of  the 
normalized  variance,  while  each  group  of  points  is  400 
seconds  of  data  on  a  given  star.  Figure  5.2  shows  very  good 
agreement  with  theory  on  this  test.  The  least-squares  fit 
has  a  slope  of  2.60,  which  is  not  a  significant  error  from 
8/3.  The  least-squa res  fit  does  not  include  the  two  data 
sets  at  high  zenith  angles  since  the  normalized  variance  is 
saturating . 

The  second  set  of  data  is  from  Albuquerque,  New  Mexico, 
acquired  when  atmospheric  conditions  were  less  than  optimal. 


Log  Secant (  Zenith  Rngle 


During  the  daylight  hours,  large  thunderstorm  cells  formed, 
dissipating  at  night  into  a  layer  of  cirrus.  About  0200, 
the  cirrus  cleared  and  we  obtained  the  data  set  in  Figure 
5.3.  The  least-squares  fit  includes  the  three  observations 
at  lower  zenith  angles.  The  slope  of  the  least-squares  fit 
is  2.48,  less  than  the  8/3  theoretical.  Based  on  the 
conditions  that  night,  the  results  are  v  .y  good. 

More  interesting  in  Figure  5.3  is  the  behavior  of  the 
normalized  variance  at  large  zenith  angles  (last  two  data 
sets) .  We  note  that  the  normalized  variance  is  no  longer 
changing  for  these  data  sets.  Clearly,  the  path  is  long 
enough  to  induce  saturation  in  the  normalized  variance.  For 
the  fourth  generation  isoplanometer  this  occurs  at 
normalized  variances  of  about  .06  to  .07.  The  10  cm 
isoplanometer  saturates  at  normalized  variances  of  .2  to  .3. 
The  aperture  averaging  has  lowered  the  normalized  variance 
between  the  two  instruments  the  predicted  amount.  Remember, 
the  calibration  constant  A'  for  che  10  cm  isoplanometer  is 
approximately  1  x  10^,  while  A*  for  the  fourth  generation 
isoplanometer  is  4  x  10  Hence,  for  a  given  isoplanatic 
angle,  the  normalized  variance  for  the  new  system  will  be  a 
factor  of  four  lower  than  the  10  cm  system.  Thus,  the 
double  annulus  system  is  saturating  at  an  equivalent 
normalized  variance.  We  still  have  not  resolved  whether  the 
increased  aperture  averaging  aids  in  forestalling 


saturation.  The  fourth  generation  system  saturates  at  its 
nominal  "saturation"  normalized  variance,  however,  it  may 
take  four  times  the  turbulence  to  make  it  saturate.  The 
fourth  generation  system  seems  to  take  reliable  data  at 
slightly  larger  zenith  angles  than  the  10  cm  system. 
Resolving  this  question  requires  simultaneous  observations 
with  both  instruments  on  the  same  stars  for  a  zenith  angle 
check.  This  test  has  not  been  possible  due  to  instrument 
availability. 

C.  APODIZED  AND  10  cm  ISOPLANOMETER  COMPARISON 
Even  though  simultaneous  zenith  angle  tests  were 

impossible  in  Albuquerque  due  to  instrument  availability  and 
atmospheric  conditions,  we  did  obtain  simultaneous  data  on 
Vega  on  August  25,  1985.  The  data  is  contained  in  Table  3. 
The  mean  isoplanatic  angle  for  the  10  cm  isoplanometer  data 
is  6.7  +  .17;  while  the  double  annulus  averaged  6.8  +  .15. 
This  shows  very  good  agreement  between  the  two  instruments 
and  removes  the  ambiguity  of  the  multiplicative  constant  in 
the  zenith  angle  test.  Table  3  also  shows  an  unstable 
night.  The  isoplanatic  angles  range  from  4.5  prad  to  9.0 
prad  over  a  one  hour  period. 

D.  ISOPLANATIC  ANGLE  MEASUREMENTS  WITH  SUPPORTING 
METEOROLOGICAL  DATA 

During  September,  1985,  the  fourth  generation  system  was 
in  use  on  Haleakalea  crater,  Maui,  HI.  One  night,  our  data 


TABLE  3 


10  cm  AND  DOUBLE  ANNULUS  COMPARISON 
Time  (Local  Standard)  10  cm  (yrad)  Double  Annulus  (urad) 


TABLE  3  (CONTINUED) 


Time  (Local  Standard)  10  cm  (_urad)  Double  Annulus  (urad) 


collection  overlapped  a  meteorological  balloon  flight.  The 
balloon  provided  altitude,  pressure,  temperature,  dewpoint 
depression,  relative  humidity,  wind  speed,  and  direction. 
Figure  5.4  shows  the  real-time  computer  output  on  17  Sep 
1985  (Star  17=Vega) .  The  solid  line  is  the  time  history  of 
the  average  intensity,  while  the  points  represent  1  second 
measurements  of  the  isoplanatic  angle.  The  spikes  in  the 
average  intensity  occurred  when  the  telescope  dome  lights 
were  turned  on.  Figure  5.4  shows  a  small  isoplanatic  angle. 
However,  no  zenith  angle  correction  has  been  applied  to  this 
data.  Figure  5.5  contains  this  15  minute  set  of  data  after 
application  of  the  zenith  angle  dependence.  The  isoplanatic 
angles  are  in  fact  fairly  large,  approximately  13  u rad . 

These  large  isoplanatic  angles  indicate  low  turbulence 
between  8-15  km.  Hence,  the  wind  profile  should  show  low 
wind  speeds  and  small  shears.  The  complete  set  of  meteoro¬ 
logical  data  is  contained  in  Appendix  B,  the  wind  shear  data 
is  as  follows.  The  balloon  data  indicates  three  shear 
layers  (2  km,  14  km,  27  km).  The  wind  speed  at  the  2  and  14 
km  levels  is  approximately  4  knots,  while  at  27  km  is  47 
knots.  These  wind  speeds  are  very  low  and  indicate  that 
very  little  shear  driven  turbulence  is  present.  Hence,  the 
large  isoplanatic  angles  agree  with  the  observed  wind 
profile  data. 
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E.  QUALIFICATION  CONCLUSIONS 

Based  on  the  zenith  angle  tests,  the  simultaneous 
measurements  with  the  10  cm  device  and  the  measurements  with 
supporting  meteorological  data,  we  conclude  that  the 
instrument  is  indeed  measuring  the  isoplanatic  angle 
correctly.  Currently,  two  fourth  generation  instruments  are 
in  operation  at  different  locations  worldwide.  This  data  is 
being  archived  at  the  Naval  Postgraduate  School  by  D.  L. 


Walters . 


As  mentioned  previously,  we  conducted  a  trip  to  Maui,  HI 
supplying  an  isoplanometer  to  the  AMOS  facility.  During 
this  trip,  weather  prohibited  acquisition  of  large  amounts 
of  data,  however,  we  did  obtain  that  data  mentioned  in 
Chapter  V,  along  with  what  I  will  present  here. 

In  this  section,  we  will  see  one  of  the  greatest  assets 
of  the  fourth  generation  isoplanometer.  That  is  its  ability 
to  make  isoplanatic  angle  measurements  effortlessly  each 
second.  The  following  plots  show  trends  in  isoplanatic 
angle  data  that  have  not  been  observed  before. 

Figures  6. 1-6.7  are  real-time  computer  output  of  each  1 
second  isoplanatic  angle  measurement.  As  in  Figure  5.4,  the 
solid  line  is  the  average  stellar  intensity  and  each  point 
is  an  isoplanatic  angle  measurement.  We  might  note  that  the 
average  intensity  is  much  more  steady  than  in  Figure  5.4. 
This  is  due  to  the  lack  of  high  thin  cirrus  clouds.  Also, 
in  Figure  6.2,  the  steps  are  caused  by  changing  the  voltage 
to  the  photomultiplier.  Note  that  these  changes  do  not 
affect  the  isoplanatic  angle  measurements.  In  Figure  6.7, 
the  deviations  in  the  average  intensity  are  due  to  the  star 
beginning  to  move  out  of  the  telescope  f ield-of-view. 

Tn  these  data  sets,  we  see  isoplanatic  angles  ranging 
from  15  Mrad  to  5  yrad.  This,  of  course,  indicates  a  very 
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unstable  night.  But,  what  is  most  interesting  are  the  low 
frequency  oscillations  present  in  the  data.  Figure  6.4  has 
oscillations  with  periods  on  the  order  of  a  one  minute, 
while  Figure  6.5  shows  oscillations  on  the  order  of  three 
minutes.  Unfortunately,  no  balloon  measurements  of  meteor¬ 
ological  data  were  performed  on  20  September.  Without 
supporting  meteorological  data,  it  is  impossible  to 
ascertain  exactly  the  mechanism  responsible  for  producing 
data  like  this. 

Based  on  the  isoplanometer  weighting  function,  it  is 

possible  to  make  some  general  comments  on  the  location  of 

the  turbulence  that  could  produce  this  type  of  oscillation. 
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Since  the  isoplanometer  is  weighted  z  ,  low  altitude 
turbulence  (boundary  layer  to  2  km)  is  not  responsible  for 
this  behavior.  Most  likely,  the  turbulence  producing  the 
oscillations  is  located  in  the  tropopause.  At  10  km,  the 
z5/3  weighting  and  Cn2  values  are  large  enough  to  produce 
the  lower  isoplanatic  angles  we  see.  Researchers  in 
Albuquerque,  NM,  indicate  that  they  observe  oscillatory 
data,  with  the  second  fourth  generation  isoplanometer,  when 
the  jet  stream  is  moving  into  the  area.  As  the  jet  stream 
approaches,  wind  speeds  in  the  tropopause  begin  to  increase. 
Periodically,  the  gradient  Richardson  number  (ratio  of 
potential  temperature  gradient  to  velocity  gradient),  will 
fall  below  the  critical  Richardson  number  (about  R.  =  1/4), 


and  the  flow  will  give  way  to  shear  driven  turbulence.  The 
jet  stream  is  not  close  enough  to  supply  the  energy 
necessary  to  maintain  the  turbulence.  Hence,  after  a  short 
period  of  time,  the  flow  becomes  stable  and  the  isoplanatic 
angles  increase.  Once  the  jet  stream  has  completely  entered 
the  area,  the  isoplanatic  angles  no  longer  oscillate,  but 
remain  low  due  to  constant  turbulence.  Thus,  one  can 
surmize  that  the  oscillatory  isoplanatic  angles  are  caused 
by  temporally  intermittent  turbulence. 

More  research  needs  to  be  done  in  the  area  of 
isoplanatic  angle  measurement  with  supporting  meteorological 
data.  Once  trends  in  isoplanatic  angle  data  have  been 
established  for  certain  atmospheric  phenomena,  balloon  data 
may  no  longer  be  necessary  for  an  accurate  picture  of 


VII.  REFRACTIVE  TURBULENCE  PROFILING 


In  the  discussion  on  the  isoplanatic  angle,  we  saw  that 

by  measuring  the  normalized  variance,  we  avoided  the  unknown 
2 

Cn  profile.  Clearly,  the  integrated  parameters  like  0g  and 

rg  do  not  provide  direct  information  on  the  distribution  of 

turbulence  with  altitude.  Many  laser  system  design  analyses 
2 

require  Cn  profiles.  However,  they  are  difficult  to  obtain 

by  remote  measurements.  In  Chapters  IV  and  V,  I  presented  a 

profile  measured  by  microthermal  probes  on  a  balloon.  This 

2 

represents  a  direct  measurement  of  C^  ,  with  very  high 

vertical  resolution  (on  the  order  of  meters) ,  along  a  single 

2 

vertical  path  through  the  atmosphere.  Actively,  Cn 
profiles  are  measured  by  acoustic  sounders  (up  to  approxi¬ 
mately  1  km)  and  pulsed  Doppler  radars  (2-30  km)  [Refs.  41, 
43],  Active  techniques  provide  good  vertical  resolution, 
but  the  cost  of  such  instruments  is  quite  high. 

In  this  chapter  and  the  next,  we  will  look  at  profiling 

2 

Cn  passively  by  direct  inversion  of  the  amplitude 
scintillation  covariance  equation  with  application  of  the 


Tikhonov  regularization  technique.  But  before  delving  into 

a  rigorous  solution,  we  should  make  note  of  the  other 

2 

passive  techniques  for  measuring  C  profiles. 


A.  PREVIOUS  WORK 


The  time  lagged  covariance  function  of  amplitude 
scintillations  for  a  point  source  and  two  equal  circular 
aperture  receivers  is  given  by  [Ref.  44], 
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where  p*  is  the  separation  between  the  detectors,  v(z')  is  a 

wind  profile,  t  is  the  time-lag,  *{K')  is  the  0.033K' 

portion  of  the  Kolmogorov  spectrum  and  D  is  the  diameter  of 

the  circular  receiving  apertures.  Often  this  equation  is 

written  for  point  receivers.  In  that  case  the  J-^xJ/x  +  1. 

2 

For  measuring  profiles,  we  let  the  time-lag,  t ,  equal 

zero,  removing  the  unknown  velocity  profile.  Hence,  B^fp') 
is  in  the  form  of  a  Fredholm  equation  of  the  first  kind, 

B  x  ( p  1  )  =  A  Jq(  z',p')Cn2(z')dz'  ,  (7.2) 

where  A  is  a  constant  and  Q(z',p')  is  the  kernel  given  by 
the  integral  over  K'.  Fredholm  equations  of  the  first  kind 
are  notorious  for  being  ill-posed  (i.e.,  errors  or  noise  in 
the  data  make  the  solution  ill-behaved).  In  fact,  Fredholm 


equations  of  the  first  kind  may  have  no  solution  at  all. 
Peskoff  [Ref.  45]  was  the  first  to  invert  this  equation  for 
point  apertures.  Shen  [Ref.  46]  realized  that  the  problem 
was  ill-posed  and  used  a  least-squares  technique  to  minimize 
error.  As  time  progressed,  other  methods  for  handling  the 
ill-posed  problem  appeared.  Heneghan  and  Ishimaru  [Ref.  47] 
proposed  a  statistical  inversion  technique  to  compensate  for 
the  large  errors  induced  by  inversion.  However,  one  should 
note  that  none  of  these  techniques  evolved  into  an  actual 
measurement  system.  People  began  to  realize  that  if  the 

2 

weighting  function,  W(z')»  (3.3),  could  be  made  to  peak,  Cn 
could  be  measured  directly  without  inversion  using  (3.1). 

1.  Remote  Probinq  with  Apodized  Apertures 


To  make  W(z')  peak  requires  filtering  in  the 

aperture  or  image  plane  [Refs.  31-33].  NOAA  [Ref.  32] 

2 

fabricated  a  Cn  profiling  instrument  (stellar  scintillo¬ 
meter)  based  upon  the  filtering  technique.  By  using  three 

2 

optical  weighting  functions,  it  measures  values  of  Cn  at 
seven  levels  between  0-20  km.  Unfortunately,  the  weighting 
functions,  W(z'),  are  somewhat  broad.  The  structure  we 
observe  in  Figure  3.2  would  not  be  as  pronounced  if  this 
profile  had  been  measured  by  the  optical  scintillometer. 
However,  the  scintillometer  does  follow  the  general  trends 


of  Cn  well. 


In  Appendix  A,  I  show  many  weighting  functions, 
W(z),  for  different  apertures  and  spectra.  This  was  an 
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effort  to  find  highly  peaked  weighting  functions  for  an 
optical  scintillometer  type  instrument.  Unfortunately,  the 
quest  proved  fruitless,  but  other  interesting  weighting 
functions  were  found. 

2.  Binary  Star  Techniques 

Several  French  researchers  [Refs.  48-50]  have 

2 

developed  a  technique  for  profiling  Cn  using  a  spacio- 
angular  correlation  function  of  two  stars  in  a  telescope 
aperture.  The  maxima  in  the  correlation  function  correspond 
to  levels  of  turbulent  layers  through  the  relationship, 


where  pi  is  the  location  of  a  correlation  maxima,  9  is  the 

angular  separation  of  the  binary  stars  and  h^  is  the 

altitude  of  the  turbulence  corresponding  to  the  maxima.  The 

2 

group  built  instrumentation  to  measure  Cn  profiles  in  this 

manner.  They  report  vertical  resolution  on  the  order  of  2 

km,  which  is  very  good.  However,  being  dependent  on  binary 

2 

stars  is  a  significant  hinderance,  esecially  when  the  Cn 
profile  vertical  resolution  depends  on  their  angular 
separation. 

B.  INVERSION  WITH  TIKHONOV  REGULARIZATION 

Jarem  [Refs.  3,  4]  implemented  a  regularization 
technique  for  handling  the  ill-posed  direct  inversion 
problem  on  the  covariance  function  for  point  detectors. 
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Jarem  defines  a  set  of  dimensionless  variables  for 
convenience , 


where  d  is  the  maximum  separation  of  the  detectors.  It  is 
interesting  to  note  that  (L/k)  is  related  to  the  Fresnel 
zone  size  since 


L,  1/2  _  , XL, 1/2  _  f 
k)  '  ^  '  (2w)1/2  ' 


where  f  is  the  Fresnel  zone  size.  With  the  new  variables, 
* (K ' )  becomes 


•  (K  • )  =  . 033K '  11/3  =  ,033(£)~11/6  K  11/6  =  (£)-11/64  (k) . 


(7.10) 


If  the  receiver  aperture  is  given  by  a  "top  hat"  function, 
I (K ' )  becomes 


I  (K  * ) 


r„,  ,k'd,  i 

2 

r  T  y/2b,i 

2J1<— 5-) 

2Ji(  2  J 

K  *D 

K1/2b 

L  ^  J 

L  2  J 

I  (K)  , 


(7.11) 


where  b  is  given  by 


b  -  D (k/L) 


1/2 


Making  these  substitutions  into  (7.8)  gives 


F  (K • )  *  2w3k2(£)-11/6*(K)I (K)  f  LdzCn2(z') 


J  LdzCn2 (z ' )cos (Kz) 


Following  Jarem,  let 


Cn2(z)  *  ir3k3(£)"11/6  Cn2(z'  ) 


(7.11b) 


giving 


VK'> 


■ 


* (K)I (K) 


/ 


dzCN  (z) 


-  / 


dzCN  (z)cos(Kz) 


2L 

k 


FX(K) 


(7.12) 


2 

(On  Cn  ,  lower  case  n  subscript  indicates  the  dimensioned 
variable,  while  upper  case  N  subscript  indicates  the 
dimensionless  variable).  Making  similar  substitutions  into 
B^tp')  leads  to 
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(7.13) 


Bx(p')  »  Bx(p) 


where 


Bx(p)  =  ^  f  dzCN2(z)  f  dKI (K ) ♦  (K)  sin2(^|) 

0 


x  J0(aK1/2p) 


(7.14) 


and 


a 


ah1'2  ■ 


Another  important  relation  is  how  the  Hankel  transform 
scales  under  the  change  of  variables.  Beginning  with  (7.4), 
using  (7.12)  and  (7.13)  gives 


2L 

k 


FX(K) 


2 if  /  Bx(p)J0  (aK1/2p)d2PdP  . 
0 


Hence,  FX(K)  is 


/  pdPBx (p)J0 (aK1/2p) 
0 
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FV(K) 


(7.15) 


With  the  relations  we  have  developed,  (7.9)  becomes 


FX(K)  *  G(K)  [2Cn2  -  Cn2(K)]  ,  (7.16) 

where 

G  (K)  =  ♦  (K )  I  (K )  . 

The  solution  to  the  integral  equation  becomes 


CN2<*>=2T  /CN2<K>e 


iKz 


dK 


(7.17) 


or 


OD 


'-F„(K> 

- — -  +  2C 

G (K) 


eiKzdK 


(7.18) 


This  is  the  ill-posed  solution  of  the  integral  equation. 
Tikhonov  proposes  using  a  low  pass  filter  function  q(K)  such 
that  the  high  spatial  frequencies  contaminated  by  the 
inversion  process  are  removed.  We  should  note  that  this 
process  does  nothing  to  errors  at  low  spatial  frequencies. 
The  regularized  solution  is  of  the  form 


(z)  = 


1 

•7- 


r 


-FX(K) 


r.  /  v  \ 


+  2C 


KM2l 


q(K)elKzdK 


(7.19) 


A  profiling  system  would  measure  B  ( p  ' ) •  Using  the 
relations  developed,  B^p')  is  scaled  to  Bx(p)  and  trans¬ 
formed  to  Fx(K).  In  (7.19),  all  the  quantities  are  knowns 
with  the  exception  of  C^2.  CN2  is  given  by 

m 

2CN2  «  /  dzCN2<z)  , 


which  can  be  defined  as  a  dimensionless  form  of  Fried's  Tq 
by 


R 


0 


t 


where 


R 


0 


4. 66ir3  .Ik  -5/6  -5/3 

sec  $  L  r0 


(7.20) 


where  *  is  the  zenith  angle.  Instruments  exist  to  measure 

rQ  which  immediately  gives  Rg- 

Another  unknown,  in  terms  of  a  practical  measurement 

system,  is  L.  In  reality,  the  value  of  L  need  only  to  be 

2 

known  approximately.  The  profile  can  be  normalized  by 

independent  measurements  of  0Q  and  rQ. 

2.  Constructing  the  Filter  Function 

Clearly,  in  any  measurement,  noise  will  be  present. 
Hence,  the  spectrum,  Fx  (K) ,  can  be  represented  by  a  "true" 


component  and  a  "noise"  component. 


9 


(7.21) 


Fx(K)  =  Fx  (K)  +  Fx  (K) 

where  the  superscript  t  and  n  imply  true  and  noise.  Hence, 

2 

CN  (K)  is  also  composed  of  a  "true"  and  "noise"  component. 

If  we  let 

Cs2(K)  =  2Cn2  -  Cn2(K)  , 

(7.15)  becomes 

FX(K)  =  G(K)Cg-2(K)  .  (7.22) 

The  kernel,  G(K),  is  type  I  to  good  approximation  in 
Tikhonov's  representation.  Type  I  implies  that  there  are 
no  zeros  along  the  real  axis  and  as  K  ♦  ®,  B  (K )  *  0. 

For  type  I  kernels,  the  filter  function,  q(K)  can  be 
represented  by 

q  (K)  =  - [G  -  ,  (7.23) 

[G(K)]Z  +  M(K) 

where  M(K)  is  a  smoothing  function.  If  the  spectral 

characteristics  of  the  noise  are  known,  and  the  spectrum 
2 

(K)  unknown,  M(K)  becomes 

M(K)  =  aK2P  ,  (7.24) 
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where  M(K)  is  a  p 


order  stabilizer  with  regularization 

2 

parameter  a.  If  the  noise  spectra  and  signal  (CN  (K)) 
spectra  are  known,  one  can  construct  an  optimal  Wiener 
filter  with 

M(K)  =  §'{ff  '  (7.25) 

where  S(K)  is  the  signal  spectra  and  N(K)  is  the  noise 
spectra . 

If  the  noise  process  is  white,  like  the  shot  noise 
in  a  measurement  system,  Tikhonov  states  that  P  =  1/2  in 
(7.24).  Following  Jarem,  we  can  make  an  approximation  to 
the  optimal  Wiener  filter,  (7.25).  Consider  an  error 
functional  of  the  form 


e  =  /  [Cn2(z)  *  CN2t(z)]2dz 


(7.26) 


2 

where  CN  (z)  has  both  signal  and  noise  components  and  t  in 

2 

the  second  term  denotes  the  "true"  component  of  CN  (z).  By 
Parseval's  theorem,  we  can  rewrite  (7.26)  as 


1 


e 


9. 


I  rc  2(k) 


2t  2 

CM ^  (K ) ]  aK  . 


(7.27) 


Now,  CN  (K)  is  given  by 


'N 


(K) 


FV(K) 

G  (K) 


[G  (K)  ]  2 

[G  (K )  ]  2  +  M  (K) 


(' 


2t 

and  CN  (K)  is  given  by  (7.16).  Substituting  (7.16),  (" 
and  (7.21)  into  (7.27)  gives 


c 


1 

17 


m 

/ 


(Fxfc(K)  +  Fxn(K)) 

cTD 


[G  (K)  ] 


[G  (K )  J  +  M  (K) 


+ 


Fxt<K> 
G  (K) 


(7 


Putting  (7.29)  over  a  common  denominator  and  performing 
appropriate  cancellations  leaves 


[- 


F  fc(K) 

G(K)Fxn(K)  -  M(K)[2CnZ  +  |-(K-y-1 


[ [G (K) ] 2  +  M (K) ] 2 


dK 

(7 


2t 

where  the  second  bracket  in  the  numerator  is  merely  CN 
Squaring  the  numerator  of  (7.30)  gives 


[G  (K)Fxn(K)  ]  2  +  2M(K)G(K)Fxn(K)CN2t(K) 
[  C  (G  (K )  ]  2  +  M  (K )  ]  2 

(M(K)C  2t  (K)  ]  2]  " 

+  - 2 - 2  dK  * 

[  [  (G  (K)  ]  ^  +  M(K)r 


Jarem  notes  that  the  middle  term, 


2M(K)G(K)Fxn(K)CN2t(K) 
C  [G  (K)  ]  2  +  M  (K)  ]  2 


=  0 


since  the  noise  spectra,  Fxn(K)  and  the  signal  CN  (K)  are 
produced  by  uncorrelated  processes.  Hence,  the  error 
becomes 


=  17  / 


[G  (K)  F  n  (K)  ]  2  +  [M(K)C  2t(K)]2 


[  [G  (K) ] 2  +  M (K) ] 2 


dK  .  (7.32) 


We  want  to  construct  M(K)  so  that  the  error  is  minimized. 
This  implies  that  the  derivative  of  the  error  functional 
with  respect  to  M(K)  vanishes  or 


dM  (K ) 


=  0 


(7.33) 


After  performing  the  di f f erent i a t ion ,  we  find  that  the 
minimizing  M(K)  is  given  by 


M(K)  = 


F  n(K) 
x  v  ' 

2t 

Cn^(K) 


(7.34) 


Since  our  noise  process  is  white,  Fx  ( K )  can  be  approximated 

2t 

by  a  constant  N.  Jarem  makes  an  assumption  about  CN  (K) 
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that  is  invalid  for  vertical  soundings.  His  data  is  along  a 

2 1 

2  km  horizontal  path.  Hence,  CN  (z)  is  basically  a 
constant  making  CN2t(K)  behave  as 

A  sin  K 
K 


Substituting  back  into  (7.34)  gives 


M  (K)  = 

A*  sinz  K 


.  2 


The  average  value  of  sin  K  is  1/2  making  M(K) 


M  (K)  =  aK  , 


(7.35) 


where 


a 


(7.36) 


This  corresponds  to  a  stabilizer  of  order  p  -  1  and 

represents  an  optimal  Wiener  filter.  It  is  interesting  to 

note  that  M (K )  is  completely  independent  of  the  aperture. 

For  vertical  paths,  CN2t(z)  varies  over  five  orders 

2t 

of  magnitude.  Hence,  the  assumption  of  a  constant  C^j  (z) 
However,  it  does  represent  a  first  order 


is  not  valid. 


approximation  and  is  a  starting  point  for  numerical 

simulation.  If  problems  occur  in  the  inversion,  one  can 
2t 

model  (z)  as  an  exponential  decreasing  in  z,  giving  a 

new  smoothing  function  M(K).  Therefore,  we  have  two 
stabilizers  of  order  p  =  1  anc  p  =  1/2  to  attempt  successful 
inversion . 

For  Jarem's  inversions,  he  finds  that  the  regulari¬ 
zation  parameter,  a,  is  nominally  10  ^  for  the  horizontal 
paths.  Chapter  VIII  investigates  this  inversion  technique 
on  the  500  m  vertical  resolution  data  of  Good  shown  in 
Figure  3.2.  But  before  moving  on  to  the  simulations,  I  will 
investigate  two  analytical  results  that  we  will  find 
convenient . 

In  (7.14),  the  z  integration  can  be  completed 
analytically.  First,  interchanging  the  K  and  z  integrations 
gives 

ao 

BX(P)  =  7  J  dKG  (K)Jq  (aK1/2P) 

0 

After  completing  the  z  integral,  we  have 

00 

=  7  /  dKG  (K) Jn  (aK1/2p)  [2C  2  -  C  2(K) ]  ,  (7.38) 


j  dzCN2  (z)  sin2  (^T")  . 

00 

(7.37) 


BV(P) 


where  the  sin  term  is  expanded  using  the  double  angle 


relation  as  before.  For  test  purposes,  it  is  convenient  to 

2  2 
input  CN  (z)  as  a  cosine  wave.  Hence,  (K)  becomes  a 

delta  function  giving 

Bx(p)  /  dKG  (K) JQ (aK1/2p)  [2Cn2  -  6(K-K')1  .  (7.39) 

0 

—  2 

If  the  cosine  wave  has  no  offset,  2CN  is  zero  and  the  K 
integral  has  an  analytic  result, 

B  ( p )  =  -  G(K')Jn(aK,1/2p)  ,  (7.40) 

X  IT  U 

where  K'  is  the  wavenumber  of  the  input  cosine  wave.  Thus, 
(7.40)  is  an  analytical  check  for  softwave. 

3.  Noise  Computations 

Detectors  in  an  instrumentation  package  that  measure 
( p ' )  will  add  shot  noise  to  the  signal.  From  [Ref.  51], 
the  average,  rms  shot  noise  current  is  given  by 

iN2  ( v )  =  2elAv  ,  (7.41) 

where  Av  is  the  frequency  bandwidth,  I  is  the  average 
current  from  the  detector  and  e  is  the  electron  charge.  I 


V  *  J  ■  -  • 


depends  on  the  intensity  of  the  starlight  and  quantum 
efficiency  of  the  detectors.  To  find  the  irradiance  of  a 
given  star  on  the  telescope  aperture,  we  can  use  the 
following  [Ref.  52]  relation, 

*1 

log  (y-)  =  •4(Mag1  -  Mag2)  ,  (7.42) 


where  1^,  I2  are  the  stellar  intensities  on  earth  and  Mag^ 

and  Mag2  are  the  magnitudes  of  the  different  stars.  Using 

the  sun  as  a  reference,  we  find  that  the  star  <*  Auriga 

—8  2 

(Capella)  has  an  irradiance  of  1.05  x  10  W/m  on  the 
telescope  for  Planck  radiation  between  0.4-0. 6  inn  (visible). 
For  the  apodized  aperture,  the  detector  sees  2.15  x  10  10  W. 
The  quantum  efficiency  is  [Ref.  53] 


Q.E. 


(7.43) 


where  ng-  is  the  number  of  photoelectrons  emitted  by  the 
detector/second  and  n^  is  the  number  of  incident  quanta/ 
second  on  the  detector.  Hence,  n  can  be  found  using  the 

q 

following  relation, 


n 


q 


JP 
h  v 


(7.44) 


125 


where  P  is  the  power/wavelength  interval  of  the  starlight  on 

the  detector,  h  is  Planck's  constant  and  v  is  the  frequency 

of  the  incident  light.  For  the  illustration  of  Capella,  ng 

is  5.4  x  10  quanta/second  given  X  =  500  nm.  If  we  consider 

a  detector  with  a  quantum  efficiency  of  10%,  I  is  8.66  x 
-12 

10  A.  Hence,  the  variance  of  the  current  due  to  shot 

-27  2 

noise  is  1.39  x  10  A  given  a  bandwidth  of  500  kHz. 

Thus,  we  can  construct  a  theoretical  signal-to-noise  ratio 
given  by 


SNR  = 


8.66  x  10 
3.72  x  10' 


=  232.63 


(7.45) 


We  can  use  this  figure  for  constructing  Gaussian  white  noise 

in  Chapter  VIII.  As  pointed  out  earlier,  the  shot  noise 
— 1/2 

increases  with  I  .  Since  the  signal  is  proportional  to 
the  intensity,  the  signal-to-noise  ratio  increases  with 

I1/2. 


4.  Inversion  Conclusions 


The  derivation  of  sections  1  and  2  is  completely 
general  in  terms  of  the  receiver  apertures.  Hence,  the 
receiver  apertures  can  be  apodized  to  optimize  inversion 
performance.  We  hope  that  the  additional  spatial  filtering 
done  by  the  aperture  averaging  will  aid  in  regularizing  the 
ill-posed  nature  of  the  inversion  process. 


VIII.  NUMERICAL  SIMULATION  OF  THE  AMPLITUDE  COVARIANCE 
FUNCTION  INVERSION  WITH  TIKHONOV  REGULARIZATION 

2 

In  this  chapter,  we  hope  to  see  if  Cn  profiling  by 
direct  inversion  of  the  scintillation  amplitude  covariance 
function  with  Tikhonov  regularization  is  a  practical 
measurement  technique.  We  can  accomplish  this  by 

2 

numerically  modeling  the  inversion  process  using  the  Cn 
balloon  data  of  Good  presented  previously.  From  this  data, 

I  construct  the  covariance  function  and  use  Tikhonov 
regularization  to  invert  the  covariance  function  regaining 
the  original  data.  What  follows  is  a  step-by-step  analysis 
of  the  numerical  simulation. 

A.  CONSTRUCTING  THE  CN2  AND  SCINTILLATION  SPECTRA 

2 

The  first  step  in  obtaining  CN  (K)  is  to  transform  Good's 

data  into  the  dimensionless  set  of  variables  described  in 

Chapter  VII.  We  accomplish  this  by  applying  (7.11b) 

2  2  2 

obtaining  CN  (z)  from  C^  ( z  * ) .  C^  (K)  is  simply  the  Fourier 
2 

transform  of  C^  (z).  For  this  test,  the  path  length  L  is 

2 

31.5  km  representing  the  first  63  points  of  the  known 

2 

profile  (500  m  to  31.5  km).  I  obtained  the  value  of  CN  at 

2 

the  surface  by  assuming  the  -4/3  dependence  of  CN  at  low 
altitudes  described  in  Chapter  IV.  Since  the  balloon  data 
has  a  vertical  resolution  of  500  m,  Az  =  0.0159.  The  spatial 


/  4 
,*  •  1 


•  ••  1 
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frequency  resolution  of  the  Fast  Fourier  Transform  (FFT)  is 
given  by  [Ref.  54]  as 


where  N  is  the  number  of  data  points.  Hence,  AK  is  0.9843. 

(We  should  remember  that  K  and  z  are  dimensionless 

2 

variables.)  The  value  of  AK  implies  that  CN  (K)  should  span 
zero  to  63.  However,  FFT's  are  based  on  the  Nyquist 
criteria,  [Ref.  55]  that  requires  two  points  per  Az 
interval.  Therefore,  half  the  frequencies  in  the  FFT  are 
duplicates,  implying  that  K  spans  zero  to  31.5. 

2 

The  analysis  in  Chapter  VII  requires  that  CN  be  an  even 

function.  This  is  convenient  since  the  Fourier  transform  of 

an  even  function  is  real.  Figure  8.1  shows  the  real  part  of 
2  2 

CN  (K).  The  modulus  of  (K)  is  simply  the  absolute  value 

of  the  real  part  since  the  imaginary  part  is  zero.  We 

obtain  the  scintillation  spectrum  FX(K)  by  applying  (7.16). 

It  is  interesting  to  note  that  the  FFT  is  very  robust  in 

the  presence  of  noise.  Using  an  algorithm  from  [Ref.  56],  I 

injected  Gaussian  white  noise  with  poorer  s igna 1-to-noise 

than  that  described  in  Chapter  VIII,  into  the  real  and 

2  2 

imaginary  parts  of  CN  (K) .  The  real  part  of  CN  (K)  with 

2 

noise  is  shown  in  Figure  8.2.  If  we  inverse  FFT  CN  (K)  with 

2 

noise,  we  obtain  the  C^  (z)  profile  in  Figure  8.3.  The 


(K)  Spectrum 
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SPATIAL  FREQUENCY  (NORMALIZED) 

2 

Figure  8.2  The  Real  Part  of  the  C  (K)  Spectrum  with  and  without  Noise 
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noise  has  very  little  affect  on  the  transform.  We  will  find 
that  very  small  amounts  of  noise  produce  drastic  results  in 
the  inversion  of  the  covariance  equation. 

B.  CONSTRUCTING  THE  COVARIANCE  FUNCTION 

I  determine  the  covariance  function  using  (7.38).  The 
numerical  integration  uses  the  Trapezoid  rule  at  each 
detector  separation  p.  We  should  note  that  the  Trapezoid 
rule  integrates  functions  with  noise  quite  well  since  it 
does  not  require  higher  order  derivatives  in  the  data. 

Figure  8.4  shows  the  covariance  function  for  a  maximum 
detector  separation  of  d  =  0.5  m.  The  detector  optics  have 
a  "top  hat"  aperture  with  a  diameter  of  1  cm.  The  form  of 
the  covariance  function  in  Figure  8.4  is  consistent  with 
theoretical  calculations  [Refs.  6,  15]. 

The  analytical  form  of  the  covariance  function  for  a 
2 

delta  function  CN  (K)  is  given  by  (7.40).  I  used  this  to 
check  the  covariance  function  software.  The  relative  error 
is  .02.  Hence,  the  covariance  function  calculations  have 
"noise"  (errors)  due  to  the  numerical  computations.  We  will 
use  this  as  a  source  of  noise  for  ill-conditioning. 

We  should  note  that  the  covariance  function  has  the 
proper  dependence  on  optics  aperture  size  and  electro¬ 
magnetic  wave  wavelength.  Figure  8.5  is  the  dependence  of 
the  covariance  function  on  aperture  size.  As  expected,  the 
covariance  function  becomes  wider  with  larger  apertures. 


NORMALIZED  COVARIANCE 


NORMALIZED  DETECTOR  SEPARATION;  ROW/D 


This  is  predicted  by  the  transform  properties  involved. 
Figure  8.6  shows  that  the  covariance  function  widens  for 
increasing  wavelength. 

9 

In  a  profiling  system  based  on  this  technique,  we 

would  measure  the  intensity  covariance  function.  Hence,  the 

covariance  function  is  the  starting  point  for  inversion  to 
2 

obtain  the  profile.  The  following  sections  describe  the 

data  manipulations  required  for  an  actual  measurement 
system.  it  is  interesting  to  note  that  experimental 
covariance  function  measurements  of  laser  beam  scintillation 
do  not  have  the  negative  portion  that  represents  the 
correlation  of  maxima  and  minima  in  the  scintillation 
pattern  [Ref.  28].  This  may  be  due  to  the  fact  that  the 
photomultipliers  measure  a  positive  definite  quantity  (i.e., 
intensity) . 

C.  CONSTRUCTING  THE  SCINTILLATION  SPECTRUM  AND  APPLYING 
TIKHONOV  REGULARIZATION 

The  scintillation  spectrum  FX(K)  is  the  Hankel  transform 

of  the  covariance  function.  The  dimensionless  form  of  the 

2 

transform  is  given  by  (7.15).  CN  (K)  is  found  by  (7.16). 

To  compensate  for  the  ill-conditioning  present,  we  introduce 

a  low  pass  filter  of  the  form  given  in  (7.23).  This  is  the 

2 

Tikhonov  regularization.  Figure  8.7  shows  C^  (K)  after  the 

3 

regulari  zation  with  a  stabilizer  of  the  form  <*k  .  Note  the 
affect  of  the  filter  on  the  higher  spatial  frequencies. 


structed  From  the  Covariance  Function  with  Tikhonov 


Figure  8.8  shows  how  CN  (K)  changes  with  stabilizers  of 

2 

different  powers  of  K.  We  obtain  CN  (z)  by  inverse  FFT. 

2 

Figure  8.9  shows  the  true  CN  (z)  profile,  unregularized 

solution  and  regularized  solution  with  «  =  1  x  10  ^  and  a 

stabilizer.  Remember  that  a  determines  the  amount  of 

ill-conditioning  allowed  to  enter  the  solution.  The 

unregularized  solution  (dots)  bears  no  resemblence  to  the 
2 

original  CN  (z)  profile.  The  regularized  profile  with  a  =  1 

x  10  does  very  well  resolving  much  of  the  500  m 

resolution  structure  present  in  the  actual  data.  The  small 

"hump"  in  the  regularized  solution  after  the  real  data  has 

gone  to  zero  is  part  of  the  ill-conditioning.  Note  that  the 

regularized  solution  has  problems  near  zero.  It  in  fact 

becomes  negative  (not  shown)  which  is  unphysical.  This  is 

very  likely  due  to  the  frequency  truncation.  Consider  a 

"top  hat”  aperture.  We  can  Hankel  transform  the  "top  hat" 

to  obtain  the  Airy  spectrum.  Hankel  transforming  back  to 

aperture  coordinates  produces  an  oscillation  in  the  original 

aperture  function.  Figure  8.10.  This  is  due  to  frequency 

truncation.  This  same  truncation  occurs  in  our  calculation 

2 

of  the  covariance  function  and  spectrum.  Hence,  an 

2 

induced  oscillation  in  CN  (K)  transforms  as  the  large 
negative  value  in  the  space  coordinate,  z.  The  fact  that  it 
is  negative  implies  that  the  oscillation  in  the  spectra  is 
out  of  phase  with  the  other  frequencies  present. 


Figure  8.8  C.T  (K)  for  Different  Stabilizers 


Effect  of  Frequency  Truncation  on  the  "Top  Hat"  Aperture 


Lowering  the  regularization  parameter  a  allows  more  of 

the  ill-posed  nature  of  the  problem  into  the  solution. 

Figure  8.11  shows  the  true  solution  and  regularized  solu- 
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tions  with  a  =  1  x  10  and  1  x  10  for  a  K  stabilizer. 
Clearly,  the  choice  of  o  is  a  trade-off  between  structure 
and  magnitude  accuracy.  The  aperture  averaging  does  in  fact 
filter  the  higher  spatial  frequencies,  however  manipulating 
o  produces  more  favorable  results  than  changing  the  optics 
diameter . 

2 

It  is  interesting  to  note  the  affect  on  CN  (z)  for 

stabilizers  with  different  powers  of  K  and  the  same 

2 

regularization  constant,  a.  Figure  8.12  shows  CN  (z)  for 

K2,  K3  and  K4  stabilizers  all  with  a  =  1.0  x  10  In 

2 

Chapter  VII,  we  found  that  K  stabilizers  correspond  to  the 

optimal  Wiener  filter.  From  this  simulation,  I  feel  that  K3 

2 

stabilizers  offer  the  best  fit  to  the  CN  profile. 

D.  REGULARIZATION  CONCLUSIONS 

Initially,  we  were  interested  in  whether  direct 

inversion  of  the  covariance  function  with  Tikhonov 

.  2 
regularization  offered  a  practical  way  of  obtaining  CN 

profiles.  I  have  shown  that  one  can  obtain  high  vertical 

resolution  profiles  (Az  =  500  m)  with  this  technique. 

However,  covariance  function  measurements  are  very  difficult 

to  make.  The  quality  of  the  measured  covariance  functions 

will  never  be  comparable  to  those  calculated  numerically  in 


IX.  CONCLUSIONS 


The  importance  of  the  atmospheric  isoplanatic  angle  and 
refractive  turbulence  profiles  to  the  Space  Defense 
Initiative  (SDI)  cannot  be  overstated.  In  this  dissertation 
I  have  demonstrated  an  improved  isoplanatic  angle 
measurement  system  and  presented  a  technique  to  measure 
refractive  turbulence  profiles.  In  the  following  sections, 

I  will  briefly  review  each  area  and  propose  future  work. 

A.  ISOPLANATIC  ANGLE  MEASUREMENT 

We  saw  that  the  original  Walters  isoplanometer  could  be 

improved  by  apodizing  the  aperture  of  the  receiver  tele- 
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scope.  To  best  approximate  the  z  altitude  weighting 
required  for  isoplanatic  angle  measurement,  I  constructed  a 
double  annulus  aperture  function  optimized  for  a  Celestron 
C8  telescope.  The  double  annulus  aperture  weighting 
function  is  a  significant  improvement  over  that  of  the  "top 
hat"  clear  aperture  used  on  the  previous  isoplanometers . 

Sensitivity  analysis  of  the  double  annulus  weighting 
function  to  strong  low  altitude  turbulence  and  inner  scale 
changes  with  altitude  showed  no  effect  detrimental  to  a 
measurement  system.  Hence,  we  fabricated  a  fourth 
generation  isoplanometer  based  on  the  double  annulus 
aperture.  This  instrument  is  designed  for  night  use  and 
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offers  many  advantages  over  the  previous  system.  First,  the 
double  annulus  no  longer  requires  the  use  of  the  Celestron 
C14  telescope  as  did  the  "top  hat"  aperture.  This  is  a 
considerable  savings  in  cost  and  instrument  managability. 
This  instrument  also  has  the  ability  to  make  isoplanatic 
angle  measurements  once  per  second  giving  data  with  very 
high  temporal  resolution.  Trends  in  the  isoplanatic  angle 
are  emerging  now  that  have  never  been  seen  before  and 
indicate  that  this  instrument  may  be  an  extremely  effective 
tool  to  probe  the  tropopause. 

In  Chapter  V,  I  detailed  tests  proving  that  the 
instrument  is  indeed  measuring  the  isoplanatic  angle 
correctly.  And  lastly,  I  presented  data  from  Maui  showing 
interesting  trends  in  the  isoplanatic  angle.  Currently,  two 
fourth  generation  isoplanometers  are  operating  at  different 
locations,  supporting  high  priority  programs. 

Clearly,  the  greatest  need  in  the  area  of  isoplanatic 
angle  measurement  is  a  painstaking  correlation  between  the 
structure  observed  in  isoplanatic  angle  measurements  and 
meteorological  data.  One  hopes  that  the  isoplanometer  might 
effectively  measure  turbulence  at  the  tropopause.  This 
would  save  money  in  balloon  flights,  etc. 

Also  of  interest  would  be  to  modify  the  fourth 


generation  system  to  provide  high  temporal  resolution  data 
during  the  day.  This  requires  the  use  of  high  quality 
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telescope  mounts  and  drives,  but  I  feel  daytime  data  is  well 
worth  the  expense. 

B.  REFRACTIVE  TURBULENCE  PROFILING 

In  this  section,  we  investigated  the  inversion  of  the 
scintillation  amplitude  covariance  equation.  I  included  the 
filtering  done  by  the  receiver  optics  and  obtained  a  closed 
form  solution  to  the  inversion  problem.  Unfortunately, 
inverting  the  covariance  function  is  ill-posed  in  a 
mathematical  sense,  meaning  that  noise  in  the  data  causes  an 
unstable  solution.  Hence,  I  incorporated  a  regularization 
technique  developed  by  Tikhonov  to  compensate  for  the  ill- 
conditioning.  Tikhonov  regularization  is  a  low  pass  filter 
designed  to  eliminate  high  spatial  frequencies  most 
contaminated  by  the  ill-conditioning. 

I  modeled  the  inversion  numerically  using  an  actual 
refractive  turbulence  profile  with  500  m  vertical 
resolution.  We  saw  that  the  Tikhonov  regularization 
reproduced  much  of  the  profiles  very  well.  However,  we  note 
that  actual  covariance  function  measurements  are  very 
difficult  to  make.  Hence,  I  recommend  pursuing  the 
weighting  function  approach  to  refractive  turbulence 
profiling  described  in  Chapter  VII.  One  should  be  able  to 
obtain  suitable  weighting  functions  given  an  appropriate 
basis  set.  This  would  allow  refractive  turbulence  profiling 
by  scintillation  intensity  variance  measurements,  operating 
much  like  the  isoplanometer . 
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APPENDIX  A 


OTHER  WEIGHTING  FUNCTIONS 

Many  different  weighting  functions  can  be  generated 

using  the  technique  mentioned  in  Chapter  III.  What  follows 

is  a  brief  overview  of  the  many  weighting  functions  I 
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generated  while  trying  to  approximate  to  the  z  '  behavior 

or  find  weighting  functions  that  peak  at  a  given  z. 

First,  a  few  words  about  the  double  annulus  aperture. 

In  Figure  3.8,  we  have  noted  the  error  remaining  at  low  z. 

This  can  be  improved  by  making  the  inner  annulus  larger, 

Figure  A.l.  However,  you  sacrifice  some  performance  at 

5 /3 

larger  altitudes.  Because  of  the  z  '  weighting  for  the 
isoplanometer ,  we  felt  error  should  be  minimized  at  higher  z 
at  the  expense  of  low  z. 

I  also  tried  other  combinations  of  Bessel  functions. 

One  is  the  "wedding  cake"  aperture  function  (in  my 

nomenclature) ,  which  is  composed  of  two  "top  hat"  functions 

with  different  transmissivities,  Figure  A. 2.  Typically,  the 

2 

"wedding  cake"  aperture  gave  a  weighting  function  of  z 
below  3  km.  Figure  A. 3.  Associated  with  the  "wedding  cake" 
is  the  "three  tier"  aperture  function,  which  has  one  more 
"layer".  Weighting  function  behavior  with  the  "three  tier" 
is  basically  the  same  as  the  "wedding  cake". 


One  can  incoherently  add  spectra  of  different  aperture 
functions  since  the  correlation  distance  is  so  small  (on  the 
order  of  centimeters),  given  a  telescope  with  a  sufficiently 
large  aperture.  Using  this  technique,  I  was  able  to  make 
the  weighting  function  bend  at  a  given  z,  Figure  A. 4.  The 
weighting  function  in  Figure  A. 4  is  an  incoherent  addition 
of  the  "top  hat"  and  "wedding  cake"  aperture  functions. 
Weighting  functions  with  this  characteristic  are  also 
important . 

From  the  Bateman  manuscripts,  I  tried  some  Hankel 
transform  pairs.  If  we  consider  the  Hankel  transform  of 
f(x),  (the  spectrum),  in  the  form 

co 

f (y)  -  /  f (x)JQ(xy) (xy) 1/2dx  (1) 

0 

the  function  x-1,/2e~ax  has  the  spectrum  y1^2  (y2  +  a2)  3//2. 
The  exponential  aperture  function  has  the  weighting  function 
shown  in  A. 5.  This  weighting  function  is  also  on  the  z 
side  of  z5/3. 

Shown  in  Figures  A. 6  and  A. 7  are  the  weighting  functions 
associated  with  sine  and  cosine  spectra.  Both  have 
analytical  aperture  functions  that  can  be  found  in  the 
Bateman  Manuscript  [Ref.  34]. 

This  is  by  no  means  a  complete  list  of  all  the  apertures 


I  found  weighting  functions  for.  However,  these  are  the 
more  interesting  and  applicable. 


WEDDING  CAKE  APERTURE  FUNCTION 


Figure  A.l  Wedding  Cake  Aperture  Function 


Figure  A.  3  Error  from  z  ' J  for  a  Wedding  Cake  Aperture  Function 


APPENDIX  B 


MAUI  METEOROLOGICAL  BALLOON  DATA 

LO-CATE  WO-8000  UPPER  AIR  DATA  MANAGEMENT  SYSTEM 
REV2.11  (C)  1984 

DATE (GMT)  TIME (GMT)  ASCENT  NO.  FLIGHT  EQ.  STATION 


850917 

04: 16 

17 

5 

0 

TIME 

PRESSURE 

HEIGHT 

TEMP 

DP-DEP 

MIN. 

MBARS 

M-MSL 

CEL. 

CEL. 

0.00 

711.  4 

0 

12.  0 

21.  1 

0.65 

700.0 

135 

11.  2 

17.8 

1.73 

682.  4 

348 

11.  2 

14.  6 

2.73 

659.4 

633 

9.5 

15.0 

4.  00 

627.7 

1038 

4.6 

12.  0 

6.  10 

577.7 

1709 

0.0 

11.2 

6.  15 

576.4 

1727 

-0.  1 

11.  3 

6.  28 

573.2 

1788 

90.8 

17.8 

6.  35 

571.7 

1816 

-0.  3 

11.  1 

7.80 

540.6 

2261 

-4.1 

7.1 

8.75 

518.8 

2585 

-5.  5 

12.  6 

8.82 

517.  3 

2608 

-2.7 

12.8 

8.  88 

515.  7 

2632 

-5.  9 

12.  3 

9.67 

500.0 

2874 

-7.  3 

13.1 

9.88 

495.  7 

2941 

-7.  9 

15.0 

11.40 

466.0 

3418 

-11.6 

11.5 

12.  80 

443.  3 

-12.  6 

14.  5 

15.55 

400.0 

4576 

-18.  4 

14.0 

22.  23 

305.  6 

6530 

-32.  5 

10.  3 

22.65 

300.0 

6660 

-33.  5 

10.0 

P106277 

SONDE  NO. 
6685 


RH 

PERCENT 

22.0 
28.  3 
36.  2 

34.7 

41.7 
42.9** 

42.  6 
46.9 

43.  3 

57.7 

36.8 
37.1 

37.6 

34.6 
29.0 
38.4 
28.  9 
28.3 
35.  1 

35.6 


26.  60 
27.87 
0 

31.30 
34.  78 

36.  72 

37.  13 
37.40 

38.  40 
38.45 
38.  58 
40.38 
43.  53 
44.83 

5.  3 


45.  63 


56.  53 


59.  42 
63.  05 
66.  17 

66.75 


70.  00 
70.07 
70.13 


PRESSURE  HEIGHT 

MBARS  M-MSL 

250.0 

7913 

233.8 

8359 

216.  1 

8876 

200.0 

9377 

169. 1 

10427 

153.2 

11022 

150.0 

11147 

148.0 

11227 

140.  6 

11531 

140.1 

11554 

139.0 

11606 

126.9 

12143 

109.  0 

13029 

102.  3 

13396 

100.  0 

13527 

98.4 

13620 

87.4 

14310 

70.0 

15633 

64.9 

16092 

57.  5 

16830 

51.  6 

17495 

50.0 

17690 

41.  6 

18835 

30.0 
20.  0 
19.6 


20883 

23427 

23554 

27711 

28794 

28794 
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28794 


TIME 

MIN. 

HEIGHT 

M-AS 

DIRECTION 

DEGREES 

SPEED 

KNOTS 

HEIGHT 

FT-MSL 

20 

5877 

234 

10.  6 

19281. 6 

21 

6169 

257 

13.  1 

20240. 8 

22 

6462 

256 

19.3 

21200.0 

23 

6770 

248 

21. 1 

22210. 5 

24 

7083 

246 

26.9 

23239.4 

25 

7397 

265 

21. 1 

24268. 4 

26 

7715 

284 

25.6 

25311. 6 

27 

8054 

292 

33.  5 

26423. 3 

28 

8395 

295 

37.6 

27541. 5 

29 

8662 

296 

39.0 

28418.8 

30 

8943 

291 

35.3 

29339.8 

31 

9277 

283 

31.0 

30435.6 

32 

9588 

268 

25.6 

31456.6 

33 

9889 

255 

26.  2 

32445. 6 

34 

10191 

240 

27.2 

33434.6 

35 

10494 

239 

28.  0 

34428.0 

36 

10801 

229 

37.7 

35437. 7 

37 

11107 

246 

36.  0 

36440. 2 

38 

11409 

260 

29.0 

37432. 3 

39 

11730 

284 

20.  1 

38485. 2 

40 

12029 

247 

23.4 

39464.0 

41 

12316 

243 

21.9 

40408. 2 

42 

12598 

248 

18.4 

41331.0 

43 

12879 

255 

16.  7 

42253.  3 

44 

13161 

251 

13.0 

43178.  2 

45 

13443 

209 

11.8 

44103.  5 

46 

13723 

197 

14.4 

45021.5 

47 

14002 

203 

8.  3 

45939. 2 

48 

14282 

24 

2.9 

46856. 9 

49 

14575 

50 

10.  9 

47816. 8 

50 

14868 

79 

15.7 

48781.4 

TIME 

MIN. 


HEIGHT 

M-AS 


DIRECTION 

DEGREES 


SPEED 

KNOTS 


HEIGHT 

FT-MSL 
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